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FOREWORD 

• 

It gives me considerable pleasure to write a foreword 
to this book by a colleague who has given me much 
assistance, as we have together made a mathematical 
attack on some rather difficult engineering problems. 

Progressive development in electrical engineering 
has been rapid, partly because many of the problems 
encountered have been such as to yield to mathe- 
matical analysis. This has made accurate prediction 
of the behaviour^ of electric circuits possible. As 
electrical engineering develops, circuits become more 
complicated and the growing importance of transient 
effects makes problems harder to solve. It is evident, 
therefore, that the electrical engineer must learn to 
use more powerful mathematical tools. One of the 
most useful for the analysis of electric circuits is thp 
Heaviside Operational Calculus. This is, in effect, 
a simple shorthand technique for evaluating the results 
of Fourier integral analysis. 

This Monograph was prepared at the suggestion of 
Professor Willis Jacksoft, D.Sc., D.Phil., M.I.E.E., 
of Manchester University, with the needs of electrical 
engineering students in mind. It is based on an 
‘ Out-of-Hours * course of lectures delivered by Mr. 
Josephs to engineers at the Post Office Research 
Station. The scheme adopted in the original lectures 
and closely followed in this book is to base electric 
circuit theory on a far-reaching theorem which, so far, 
appears to have escaped the notice of engineers. This 
theorem has been re-constructed from the scattered 
papers of Heaviside and is probably the last he ever 
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deduced ; for this reason the author has called it 
‘ Heaviside’s Last Theorem \ The famous Carson 
integral equation may t be obtained as a corollary to 
this theorem. It also leads to many others involving 
Fourier and Bessel integrals, elliptic functions, &c., 
which have a wider field of application than electric 
circuits. 

It is hoped that this little book will be of service 
in preparing the ground for a greater use of Heaviside’s 
Operational Calculus and that it will lead to a fuller 
and more just appreciation of the great work done by 
Heaviside in connexion with the foundations of electric 
circuit theory. 

G. RADLEY 

POST OFFICE ENGINEERING 

RESEARCH STATION 
October 1^44 
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CHAPTER 1 


THE FOUNDATIONS OF ELECTRIC CIRCUIT 
THEORY 

INTRODUCTION 

The object of making an analysis of an electric circuit 
problem, stated in a most general way, is to formulate as 
a function of the time the current which flows in any part 
of an electric network in response to a suddenly applied 
voltage of arbitrary form. This means that the analysis 
is not limited to a discussion of steady-state oscillations 
only, but includes the study of transient phenomena 
following abrupt circuit changes. Whilst a knowledge of 
the symbolic steady-state analysis, involving the operator;, 
is sufficient for the solution of many important technical 
problems, it is quite inadequate for the solution of many 
problems which now face the electrical engineer. In tele- 
graphy, for example, it is well known that the signalling 
speed is always limited by the transient effects ; while in 
telephony, the transient effects form the limiting factor to 
the distance over which telephonic communication can be 
established. Consequently the analysis and design of 
electric networks with special reference to their behaviour 
in the transient state and when subjected to overlapping 
transient forces are of considerable importance. There 
can be little doubt that the best method of attacking the 
more difficult problems of electric circuit theory is the 
direct operational method invented by the late Mr. Oliver 
Heaviside. In this book an endeavour has been made to 
review and discuss in simple language the application of 
Heaviside’s operational calculus to electric circuit problems. 

CIRCUIT PARAMETERS 

When tackling a circuit problem for the first time, 
Heaviside began by writing down the descriptive differ- 
. t 
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ential equations of the problem, which he could readily 
obtain from known physical laws. The next step lay in 
the solution of these equations, and it was here that 
Heaviside found the ordinary classic methods of solving 
differential equations presented formidable mathematical 
difficulties and in certain important cases broke down 
completely. This was the reason why he parted with 
academic methods and invented mathematical procedure 
of a revolutionary kind. 

We begin by discussing the mathematical basis of circuit 
theory. A circuit may be defined as a physical system in 
which the varying quantities may be expressed in terms 
of time and some other dimension. The characteristics of 
such a system may be expressed in terms of circuit para- 
meters. In the electric circuit there are three classes of 
parameters with which we have to deal. These parameters 
may be defined as follows : 

(1) Definition of the Dissipative Parameter 

The dissipative parameter of an electric circuit is called 
the resistance. Resistance is defined by the equation, 
e - Ri, which states that the voltage e across a pure resist- 
ance is proportional to the current i through it ; the 
proportionality factor being the value of the resistance R. 
If e is in volts and i % in amperes, the resistance R will be 
in ohms. 

(2) Definition of the Inertia Parameter 

The inertia parameter of an electric circuit is called the 
inductance. Inductance is defined by the equation, 

e » Ljg, which states that the voltage e across a pure 

inductance is proportional to the rate of change of current 
through it ; the proportionality factor being the value of 
the inductance L . Thus a coil has unit inductance if unit 
rate of change of current through it produces a unit voltage 
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across it. If e is in volts, i in amperes and t in seconds, 
then the inductance L will be in henrys. 


(3) Definition of the Spring Parameter 
The spring parameter of an electric circuit is called the 
capacitance. Capacitance is defined by the equation, 
de 

i - Cj-, which states that the current i through a pure 

condenser is proportional to the rate of change of voltage 
across it ; the proportionality factor being the value of the 
capacity C. Thus a condenser has unit capacity if unit 
rate of change of voltage across it produces a unit current 
through it.. The above definition may be written in the 
form 



which states that the capacity C is given by the charge q 
required to produce unit voltage across it. If e is in 
volts and the charge q in coulombs, the capacity C will 
be in farads.* 

In order to fix ideas, consider a simple electric circuit 
consisting of a resistance i?, an inductance L and a capacity 
C, in series. Let i denote the current in the circuit and 
e the applied voltage. Then the potential drop across the 

di 

resistance is iR, the inductance drop is L if and the drop • 


across the condenser- is g, where q is the charge on the 
condenser. 

Now from Kirchhoffs law relating to the potential 
drop around the circuit, we may write 


iR+ 4 +q c 


e 


* For a discussion of dimensions, see Physical Constants by 
W. H. J. Childs (Methuen’s Monographs on Physical Subjects). 



multiplying both sides of this equation by i, we get 

. tD d u 2 d f . . 
t 2 R+ jT. - ~+ « 

dt * 2 2C 


since i - The first term M expresses the rate at 

which electrical energy is being converted into heat. The 

second term ^ ^ expresses the rate of increase of the 

0 2 

magnetic energy. The third term expresses the rate 


of increase of the electric energy. The right-hand side 
is clearly the rate at which the impressed force is delivering 
energy to the circuit. Thus we may regard a resistance 
element as a device for converting electrical energy into 
heat. An inductance element may be regarded as a device 
for storing energy in the magnetic field. A capacitance 
element may be regarded as a device for storing energy 
in the electric field. 


kirchhoff’s laws 

There are two distinct problems involved in circuit 
analysis. The first is the formulation of the descriptive 
differential equations from known physical laws, and the 
second is the mathematical solution of these equations. 
The descriptive equations of a circuit problem may be 
established in a number of different ways. For example, 
they may be based on Maxwell’s dynamical theory. There 
can be little doubt, however, that the simplest basis for 
the equations of electric circuit theory are to be found in 
Kirchhoff’s Laws. For the electric circuit these laws 
state that: 

(1) In any closed circuit in a network of conductors the 

algebraic sum of the potential differences is zero. 

(2) The algebraic sum of the currents which meet at 

any point in a network of conductors is zero. 

The law which states that the sum of all the voltages 
about a ‘closed circuit Is zero follows from the principle 
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of the conservation of energy. All voltages must be 
included, not* only the applied voltages, and the voltages 
due to the circuit elements, but* any voltages induced by 
coupling with neighbouring circuits. The law which states 
that the sum of the currents flowing into a junction point 
is zero follows from the principle of the conservation of 
electricity. 

As an example of the appropriate mode of setting up 
the circuit equations, consider the two-mesh network 
shown in Fig. i. 



Let e(t) represent the applied voltage as a function of 
the time t. Let the two-mesh currents be represented 
by i x and i % respectively. Then, if the positive direction 
of the current in each mesh is taken to be clockwise, the 
resulting current in the mutual branch (consisting of a 
resistance R xl , an inductance L n and a capacity C 12 ) can 
be written down from KirchhofFs second law. Applying 
Kirchhoff’s first law to the two meshes, we obtain 

e(t) - |(*„ + R n )h + (£„ + <*} 

- + LxTTfo + C^j’l dt ) 
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O » -|(i?22 + Rn)h + C^z2 + + + 

~{ R ^ i+Lu it + u;,\ iidt } 

K we let p denote the operator ± and p ~ 1 the operation 

j dt } the above equations can be written in the Ohm’s law 
form, 

e(t) - Z n (p)i x - Z n {p)h 
o ~ - Z u (p)ii + Z n {p)h 

where the impedance functions are 

Zu(P) 10=1 (-^u + Rvt ) + p{Lu + As) + + 

Zn(p) " (i?22 + ^ 12 ) + jK-^22 + ^ 12 ) + P + 

% n ( P ) “ ^12 + 

It will be observed that our definition of a circuit includes 
mechanical systems not usually regarded as circuits. On 
account of the similarity of their descriptive differential 
equations, however, they may be treated as such. All 
rigid mechanical oscillating systems of a single degree of 
freedom may be treated as mechanical circuits, together 
with coupled combinations of such systems. This includes 
such systems as oscillograph vibrators, telephone dia- 
phragms, rotating shafts and deflecting beams, &c. 

the canonical equations 

We proceed now to a discussion of the canonical equa- 
tions of electric circuit theory. Our reasons for doing this 
are to show how the ordinary method of solving sets of 
simultaneous differential equations breaks down in the 
general case ; and to show how Heaviside overcame the 



* ELECTRIC CIRCUIT THEORY 7 

fonnidable mathematical difficulties associated with classic 
analysis by attacking the problem from a new viewpoint. 

An electric network is a connected set of circuits or 
meshes each of which can be regarded as made up of 
resistance, inductance and capacitance elements. Suppose 
the network has n independent meshes numbered from 
i to », and suppose the corresponding mesh currents be 
denoted by i Xt h • • • up to i n . Let R gg , L gg) C gg represent 
the total resistance, inductance and capacity in series in 
mesh £, and let R gk , L gh , C gh represent the mutual resist- 
ance, inductance and capacity between the circuits g and h. 
Kirchhoff ’s first law enables us to write the circuit equations 
that must hold for each and every one of the n meshes. 
Thus, if e(t) represents the voltage applied to the terminals 
of the first mesh, we may write 

e(t) - i?n*i + L n -^ + dt = Z xx (p)i x 

but this circuit is coupled to a second circuit in all three 
ways, so we may write 

+ ^)‘ 

e[t) “ Z n (p)i x + Z Xi (p)i t 

If there are n coupled circuits in the network, there will 
be n such terms, and we may write 

€ (t) = Z xx (p)i x + Z u (p) h + . . . + ^i(n — 1 )(p)i n — 1 + Z\n(P)in 
If Z n (p) represents the impedance function of the second 
mesh, then we may write for the second mesh 

o - Z n (p)i x + Z n (p)i s +*>..+ Z 2 ( n -i)(p)i n ..i + Z 2n (p)i n 
A similar equation may be written down for the other 
meshes of the network. Hence the complete equations 
for the general case may be Written in the form 

e (t) « Z n (p)i x Z ln (p)i n 

o - Z nl (p)h + . . . + Z nn (p)i n 



e(t) - Z n (p)i x 


^12 


+ PL 18 
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These equations, together with a description of the initial 
condition of the network, completely specify its perform- 
ance. Their solution agd interpretation form the subject 
of electric circuit theory and in the study of their solution 
we shall find the most direct introduction to Heaviside’s 
operational calculus. It should be noted that the canonical 
equations, (i) are capable of immediate generalization. 
The variable to be determined need not be the current. 
It may equally well be any other linear variable with which 
we may happen to be concerned. It must be noted that 
the circuit parameters are fixed and that R gh = R hg , 
L gh “ L Iig and C gh » C hg . This reciprocal relationship 
means that there are no concealed sources or sinks of 
energy in the network. To solve the set of equations (i) 
means that we must deduce mathematical expressions for 
the currents which substituted into the right-hand members 
will give the left-hand members. These expressions must 
also satisfy the terminal conditions at reference time t - o. 

THE STEADY-STATE SOLUTION 
In taking up the mathematical solution of the canonical 
equations (i) we shall start with the exponential solution. 
This is of fundamental importance to the engineer because 
it serves as the basis of the steady-state solution of alternat- 
ing current problems. Thus we consider exponential 
currents and voltages of the form exp {jcot), where jco has 
its usual significance in alternating current theory. If we 
assume that the applied voltage e(t) * we know that 
the expression for the resultant current in each branch of 
the network will contain an exponential steady-state term 
of the same form. This follow's because an exponential 
does not change in form when differentiated or integrated. 
Consequently the set of differential equations (i) may be 

reduced to algebraic equations because g (or p) may be 
replaced by jw and f dt (or p~ l ) replaced by (/&>)“*. 
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Having 1 algebrized ’ the differential equations (i) it is 
now a simple matter to determine any steady-state variable 
Ik by the elementary methods applicable to first-degree 
algebraic equations. Hence 


T M 1K {j(o) 

K= D{j<o) 


(*- I, 2, ... w) 


where M 1K (jao) is the minor determinant of the first row 
and Ktb. column, and 


'D(jco). 


Z n (j(u) . . . Z ln (jco) 
Z nl (jco) . . . Z H „(joj) 


is the determinant of the coefficients. 

It is an easy matter to put 1 K in an Ohm’s law form, thus : 


Ik-i/Z&o) 

where the impedance function Z(jo>) is the ratio of the two 
determinants D(jcu)/M lK {j(o ). 

Now assume that the applied voltage varies as E cos cot. 
Writing E cos cot in its exponential form we have 

j £ 

E cos cot m — + e~^ wt ) 


Consequently the applied voltage may be considered in 
two parts, each of which can be applied separately, and 
the results added, because we are dealing with a linear 
network with fixed parameters. Using the Ohm’s law 
form for I K we have 


T Eel* 1 Ee-W 

/K “iz0aT) + ^FF) 

In general, Z(jco) and Z{-joo) are conjugate complex 
numbers, and consequently 

Z(±jco)-\Z(jco)\e±l+ 

Hence we obtain 


2 
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which can be written as 

This in concentrated form contains the whole theory of 
the symbolic solution of alternating current problems. It 
will be observed that only steady-state currents are 
involved and such currents are ‘ forced ’ currents. They 
vary with the time in precisely the same manner as do 
the electromotive forces. Such currents are, however, 
only part of the total currents ; the complete solution must 
contain transient as well as steady-states. In mathe- 
matical language the complete solution of the canonical 
equations (i) must include both ‘ particular ’ and * comple- 
mentary ’ solutions. 

If it is found that the steady-state solution of a problem 
is sufficient, then the need for the use of the Heaviside 
calculus does not arise. In cases such as this the well- 
known symbolic method outlined above is sufficient. 
Perhaps it is now generally recognized that it is not possible 
to predict completely the performance of an electric network 
from a single steady-state solution. The growing practical 
importance of non-periodic phenomena demands a know- 
ledge not only of the steady-state but also of the transient 
response. Both of these result simultaneously from an 
operational solution. 


THE TRANSIENT SOLUTION 

In order to obtain the transient terms we must find the 
complementary function which satisfies the canonical 
equations (i). Let us assume that a solution of .these 
equations exists in the form * 

h m F n e* h - F tl e H , . . . &c. 


The correctness of this assumption may be checked by 
substituting into the set of equations (i) and putting 
e(t) equal to aero. If we make these substitutions and 
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then divide out the exponential the following set of 
algebraic equations is obtained : 


o m Z u (X)F n + . . 

• + Z ln (X)F ln | 

o~Z nl (X)F nl + . 

. . + Z nn (X)F n J 


Such a set of algebraic equations has a solution only when 
the determinant of the system vanishes ; that is, when 
D(X) = o. The equation D(X) « o is known as the deter- 
minantal equation of the network ; and if X x> A a , A 3 . . . &c., 
are the roots, then transient expressions of the form 
i g » F gl e ?it + F y2 e ?it + ... See. 
will satisfy the canonical equations (i). This is the form 
of the transient portion of the solution. The values of 
^gb Fq2i F g $ . . . See., are constants of integration which 
must satisfy imposed boundary conditions. This intro- 
duces us to the idea of terminal conditions which are of 
fundamental importajace in circuit analysis. The terminal 
conditions denote the state of the system when the e.m.f. 
is applied. The number of known terminal conditions 
required to satisfy the complementary function must be 
equal to the number of roots of the equation D(X) - o. If 
this number of independent terminal conditions is not 
known the complete solution of the canonical equations (i) 
by this method is indeterminate. 

BREAKDOWN OF THE CLASSIC METHOD 

So far we have followed the ordinary classic method of 
solving the canonical equations (i). The practical difficul- 
ties of this method of solution are : 

(i) The determination of the arbitrary constants Fgb 
Fg2i . . . &c., from known terminal condi- 
tions ; and 

(a) The determination of the roots X i} . .* &c., 
from the determinantal equation D(X) - o. 

With regard to (i) : in the general case it is quite im- 
possible to determine terminal conditions from the know- 
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ledge that if an e.m.f. be applied to a dissipative network 
at time t » o, then at that instant, all the currents in the 
inductances and all the oharges on the condensers are zero. 
With regard to (2) : in the general case it is quite impossible 
to determine the complex roots of the equation Z)(A) - 0. 
It was at this point that Heaviside parted with classic 
analysis and invented his operational calculus, with which 
we will now be directly concerned. 

THE HEAVISIDE METHOD 

From the preceding analysis it can be seen that Heaviside 
was concerned with the determination of the dependent 
variables i lt *„ . . . i n , of the canonical equations (1), as 
functions of the time t, with the terminal condition that 
these variables are all zero for the independent variable 
t < o. If the number of closed circuits or meshes is made 
infinite and the size of the circuit elements infinitesimal, 
the canonical equations (1) change into a partial differential 
equation of the wave-motion type. In his attack on the 
general problem Heaviside made use of the principle of 
superposition. When several electromotive forces are 
applied simultaneously to a network of fixed parameters, 
each produces its own effect independently of the others. 
Thus the effect of each alone may be calculated and the 
results added to obtain a resultant due to the simultaneous 
action of the various electromotive forces. This suggested 
to Heaviside that he could deduce a superposition theorem 
expressing the independent variables in terms of voltage 
and admittance functions which have been made to satisfy 
the necessary terminal conditions of the network. Thus 
it became apparent that these variables are capable of 

t 

expression in the form Ee(t - y>)d{h(yj)} t provided the 

0 

auxiliary variables h are dependent upon the general 
equations (1) and are all identically zero for t < o. 

Heaviside found that if he replaced e(t) in the equations (1) 
by a very simple discontinuous function, namely, a unit 
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voltage applied suddenly at reference time t - o, he could 
obtain for the prescribed boundary conditions an auxiliary 
set of equations in the auxiliary variables h l} h 2 , h 9y . . . h n . 
Thus 

[*] “ Z n (p)h x + . . . + Z ln (p)h n ) 

• ■ ( 2 ) 

0 “ Zni{P)h\ + . . . + Z nn (p)h n J 
where [i], the function on the left-hand side, written in 
accordance with the Heaviside notation as unity, is identi- 
cally zero for t < o and unity for t > o and h lt h 2 , . . . h n 
are identically zero for t < o. Thus the auxiliary equa- 
tions (2) are made to satisfy the necessary terminal condition 
of the canonical equations (1), and consequently from the 
principle of superposition it is possible to show that 
solutions for i 1} i 8 , t 3 , . . . &c., exist in the form 

t 

- v) h (y>) d v ■ • ■ • (3) 

0 

where the symbol \p represents the variable of integration 
and i(t) represents the current which will flow when the 
impressed voltage is given by the time function e{t). 
A short method of establishing equation (3) is given in 
Chapter 8. The function h(t) is called the ‘ Indiciai 
Admittance ’ and represents the current which flows in 
the mesh considered when a unit voltage [1] is applied to 
the input terminals of the network. It can be seen that 
the superposition theorem expressed by equation (3) 
mathematically relates the current in any mesh to the type 
of applied voltage and to the parameters and connexions 
of the system. Thus Heaviside showed that the solution 
of the canonical equations (1) depends entirely on the set of 
auxiliary equations (2) and the superposition theorem (3) 

THE OPERATIONAL EQUATION 

The important problem of determining h - h(t) remains, 
and Heaviside’s method of attack was to reduce the 
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differential equations (i) formally to algebraic equations 
by replacing the differential operator d/dt by the symbol 

p and the operation J dt by p" 1 . He also put the applied 

voltage e(t) equal to [i] and thus obtained auxiliary equa- 
tions from which he could readily obtain a symbolic 
solution for a particular variable. This symbolic solution 
he called the ‘ operational equation ’ of the problem. 
From the foregoing it can be seen that the operational 
equation of a problem is the full equivalent of the descrip- 
tive differential equations, and must, therefore, contain the 
information necessary to the solution provided the signi- 
ficance of the symbolic operator p can be correctly deter- 
mined. The operational equation of a problem obtained 
4 >y algebraic processes from equations (2) may be written 
formally as 

( 4 ) 

(K> i, 2, 3, . . . n) 

where h(t)> the indicial admittance, represents the variable 
to be determined and Z(p) represents a generating imped- 
ance function containing the p operators. 

The unit function [1] represents a discontinuous function 
of the time which is zero until t equals zero and unity 
thereafter. This unit function is impor ant because it is 
the simplest possible discontinuous function and is always 
tacitly assumed in electrical theory. As an example of 
this fact, suppose it is required to find the current response 
i of a circuit containing a resistance R and an inductance 
L in series, to which an alternating voltage E cos cot is 
applied at reference time t - 0. The differential equation is 

+iR-E cos cot 

for which we obtain by the ordinary method 
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Rt 

ERe L 




i « sin (cot + <!>)- 

Vm + w'W R 1 + ^L 1 


where 


( f > » tan - 


±1*. 

<oL 


Now suppose that E - ioo volts, R - i ohm, L - o-oi henry 
and <o - 1 007 r, then the expression for i shows that 
o*oo2 seconds after the switch is closed the current is 
17*063 amps. The expression also shows that 0*002 seconds 
before the switch is closed (i.e. at the time of - 0*002 seconds) 
the current is - 20*9 1 2 amps. This absurd result is avoided 
by making an assumption that the time t must be positive. 
Yet there is nothing in the mathematics from the setting 
up of the differential equation to its solution to justify this 
assumption. 

The fact is that the differential equation has not been set 



Blf continuous Cosine ■ Continuous Cosme 


Fig. a. — Continuous and discontinuous cosines, 
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up correctly ; the applied voltage is not ioo cos cot but 
ioo cos co*[i], where [i] is the unit discontinuous function ; 
the effect of operating on ioo cos out by [i] is illustrated 
in Fig. 2 . If the differential equation is set up correctly 
the factor [i] will then appear in the solution for i, thus 
showing mathematically that i is zero until t is positive. 
However, all trouble is avoided by assuming that t is always 
positive, and this simple example is sufficient to show that 
the existence of a unit discontinuous function [i] is always 
tacitly assumed in alternating current theory although not 
always explicitly stated. Since the unit function [i] occurs 
on both sides of the operational equation ( 4 ) its presence 
need not be explicitly stated and in general it will be 
omitted from operational equations. 



CHAPTER 2 


THE EXPANSION THEOREM 

AN EXPERIMENTAL DERIVATION 

In order to interpret his symbolic or operational equa- 
tion (4), Heaviside at first proceeded in an experimental 
manner. He directly compared the operational equations 
for certain specific problems with explicit solutions 
obtained by classic methods. Thus he was led to assign 
a definite meaning to his operators and he obtained by 
induction rules for solving the operational equations for 
those problems in which the classic methods failed him. 
A little later, however, he discovered certain relationships, 
which enabled him to proceed deductively as well as 
inductively. One such relationship was his famous Expan- 
sion Theorem. Heaviside established this theorem in a 
number of different ways, one of which will be reviewed 
here. 

Consider the operator i/Z^p) of equation (4) ; from 
the set of auxiliary equations (2), we may write 

I M 1 K {p) 

ZjAP)~ D(p) 

where Mi K (p ), the minor of the first row and iTth column, 
is of lower order than D(p) ) the determinant of the auxiliary 
equations (2). 

Let Ai, A a , . . . A n , denote the n roots of D(p) - 0. 

- 

Now assume that in general, 

m + , An 

m m p-h p-h p-b 

17 
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where A h A lt . . . A w are constants and independent of p. 

M(p) « A r (p- Wp - A*) . . . (p - A n ) + 
A*{p-\)(p-l*) . • . (p - A n ) + 


“ ^i)(P - - A n _j). 

Since both sides of this equation are identical for all values 
of p we may substitute l x for p throughout ; this gives 

M{h) « A^h - A 2 )(A X - A 3 ) . . . &-AJ 

Differentiating D(/>) with respect to />, we get 

(/> “ ^i)(p - A 3 ) . . . (p - A n ) + 


(P “ Ai )( P ~ A a ) • . • (/►-An-i) 

Hence D (p) ® (Aj - A 2 )(A! - A 3 ) . . . (Aj - A n ) 
Thus we obtain 

In like manner we obtain 

A t ~ M(X 2 )/D'(X t ) 


A^M^/DXK) 

Thus, when all the roots are unequal, we have 

M{P) Wi) i MOn) i 

m D'(kYp ~ A, D'(A,)-FT, + • * • + 
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So that 

i K y m<M 

z iP) - A^)D'(A k ) 

which is an explicit operational expression for the indicial 
admittance function h K (t). This expansion can be ex- 
pressed in a more convenient form. For put, as in 
Chapter i, 

Z{X)-D(1)/M{\) 

0 dZ(X) M(X)D'(X) - D(X)M'(k) 

b0 dl m [M(l)Y 

Now 2 )(A X ) “ 0 

Z\X k )-D\X k )/M(X k ) 

Thus we may write 

K 

- hK{t) ~E(p-^m K ) 

It will be observed that each term of the above expansion 
yields a linear differential equation which may be written as 


i i 

(FXj 


(r= I, 2, 3 ... n) 



Expanding the last term on the right-hand side of this 
equation by the binomial theorem, we get 


“r- 


ZWf 


Vtt r) 




The question now arises : What significance is to be given 
to the operator acting on the unit operand [i] ? 
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Assuming that the usual scheme of integration holds for 
this operator, we have by definition 


-> - j dt 


r r 2 


p-'.p - 2 — 

2-3 

from which it can be seen that 

- t n 

p ' n ~T\ ■ 

when n is some real number. Following the natural sug- 
gestion that we substitute ?/n\ y for the operator p~ n t in 
the operational expansion for u r> we get 




•^ + 1T + -^-+ 


X r n t n 

4. _i + 

«! 


It can be seen at once that this is equivalent to 

“ , 'Xrj3 

Substituting in the expansion for h K (t) we obtain 

K~n 

m y eV 


K - 


€ A h * ^ =* 71 J 


Looking back ft. the operational expansion for i/Z(p), it 
can be teen that the last term of the above equation is 
i/Z(o) by putting p — o ; so that 

h{l) '^) + Z xMk) • • • {s) 

It will be noted that equation (5) is an explicit expression 
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for h K (t), which we arrived at experimentally by using an 
unestablished operational process, and illustrates the type 
of inductive analysis originally u§ed by Heaviside. The 
summation indicates that we are to write as many terms 
as there are roots of the equation Z(p) = o, and into each 
insert the value of one of these roots. In equation (5) it 
can be seen that the first term is obviously the steady-state 
current, and the summation is of exactly the form of the 
complementary function for the differential equations of 
the network, where the constants of integration, the cause 
of the breakdown of the classic method, are now replaced 
by explicit expressions in terms of the generating impedance 
function. Equation (5) is the famous Heaviside Expansion 
Theorem. 

TRANSIENT OSCILLATIONS IN A SERIES CIRCUIT 

As an example of the application of Heaviside’s expansion 
theorem to an electric circuit problem, suppose it is 
required to find expressions for the current i in response 
to a direct voltage E applied suddenly (at reference time 
t « 0) to a circuit consisting of a resistance JR, inductance L 
and a capacity C in series. 

From the expansion theorem we have 

y v E Jyi 

' {t) ~ z(o) + E Li K z'(X K ) 

where A 2 , . . . A n , are the n roots of the equation 
Z(p) - o, and 

ZVk)- (K- 1,2,3,... n) 

The first term E/Z( 0) is obviously the steldy-state current, 
and the summation indicates that we are to writers many 
terms as there are roots of Z(p) 0, and into each insert 
the value of one of these roots. 

The circuit equation for the series circuit is 

iR+L It +1 c\ rdt - E 
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d f 

Replacing ^ by p and \dt by p~\ we get 


where 


i(R+*pL + i/pC)~E 
i-E/Zb) 
Z(P)-R + P L + JL 


The two roots of the equation Z(p) * o are 
A x » -u.+jo) 
and A 2 = - a ~j<o 

R t „ i R* 


where 


K “ 2 £ * nd "'‘IC 4Z, 2 


Differentiating Z(A) with respect to A, we get 


9 Z(A) 

dx~* 


L- 


A 2 C 


Introducing the values of A x and A 2 we get 


aZ(A)~ 

dk 


jjcoL 


ja-j, C(-x+jcoy -a +jco 
Similarly by replacing co by - co, we get 


az(zy 

3 A 


■ L - 


- zjo)L 


a -a, C(- a -;w) 2 -a-jw 


Hence A 2 Z (Ai) = 2jo)L 

snd AjZ (A 2 ) “ ~ 2 ]Q)L 

For A » 0 we have (i/AC) - 00, and consequently Z(o) - 00. 
Thus the first term E/Z(o) of the expansion theorem 
becomes E/ 00, showing that the steady-state current is 
zero. Substituting the above results in the expansion 
theorem* we obtain 


gjul e -jtd' 
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It can be seen that this equation represents. a damped 
oscillation, where the damping factor a - R/2L, and the 
frequency is given by 


= -i [± 

W L( 


jf _ > 

LC 4Z, 2 


1 R z 

If, however, j-g < ^- 2 , we notice that co now becomes 

a pure imaginary. Consequently, we put co»;co lt in the 
expression for i(t), and we obtain 




sin jojit sinh ouj 


Hence i - — ~ e ~ at sinh a> x t 

1 R 2 

Thus it can be seen that when j-g < — p, the current 

diminishes asymptotically to zero, and the motion is 
aperiodic instead of oscillatory. 

1 R 2 

In the critical case, however, when » ^- 2 , the two 

roots Xi and X 2 are equal and the expansion theorem is no 
longer applicable. We may readily arrive at the expression 
for the current in this case by noting that when — > A a , 
then co —> 0. Consequently, for very small differences 
between the roots X 1 and A 2 , co is very small and sin cot - cot. 
Putting sin cot * cot in the equation for i(t), we get- 

. E ^Et 

t = ~ T e~ sU .cot = -pg-" 

coL L 


Now this equation is independent of co and consequently 
must hold also for the limiting value X 1 ** A* and co - 0. 
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The current in this case also approaches zero asymptotically. 
Collecting results we see that if 


1 jR 2 • 

LC > 4Z 5 ’ then 1 ' 


coL 


sin wt 


1 

LC 


V 2 ’ 

£ 2 


, • Et 
then -j-e~ 


l6< ~ 'hen i- 


c^L 


1 sinhci)^ 


These results are well known. The point that we wish 
to emphasize is that they were arrived at by a simple 
algebraic process without having to determine integration 
constants from known terminal conditions. 


PRACTICAL APPLICATION OF THE EXPANSION THEOREM 

It has been shown that Heaviside’s expansion theorem 
gives an explicit relationship between the impedance 
function Z(p) and the indicial admittance function h(t). 
For electric circuits, where the number of roots involved 
is small, the expansion theorem enables one to calculate 
h(t) with very little labour and with the use of simple 
mathematics. Thus the theorem has proved useful for 
analysing the transient responses of air-cored transformers 
and intervalve couplings. It can be applied to all kinds 
of coupled circuits and, in some cases, to transmission lines. 
In all these cases the roots of the determinantal equation 
are easy to locate. When these roots are not easy to locate, 
that is, when the network has four or more distinct circuits, 
we run into trouble. This trouble is not due to any fault 
of the theorem itself, but simply because the engineer 
encounters practical difficulties when finding the roots of 
higher-degree algebraic equations, especially when some or 
all of the roots are complex. Thus great labour is involved 
in applying the expansion theorem to lumped networks of 
more than four degrees of freedom. 

It must be noted that the partial fraction expansion upon 
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which the expansion theorem solution depends imposes 
certain conditions upon the impedance function Z(p). 
Thus the impedance function must have no zero root and 
no repeated roots. In all passive electric networks these 
conditions are satisfied. The case of equal roots, which 
may occur when the network’ contains an internal source of 
energy, such as an amplifier, can be treated by assuming 
unequal roots and then letting these roots approach equality 
as a limit. In general, the expansion theorem gives the 
indicial admittance function (whenever this function is 
capable of such expression) in terms of the normal or 
characteristic vibrations of the circuits concerned. 


3 



CHAPTER 3 

EXTENSION OF THE EXPANSION THEOREM 


THE INDICIAL ADMITTANCE FUNCTION 


It was stated in Chapter i that by virtue of the super- 
position theorem expressed by equation (3), the problem 
of circuit analysis was reduced to the determination of the 
indicial admittance function h(t). In other words, the 
mathematical formulation of the current flowing in any 
part of the network in response to a unit voltage (zero 
before, unity after time t - 0) is of fundamental importance. 
It was shown in Chapter 2 that a formal expression for the 
indicial admittance function is furnished by the Heaviside 
expansion theorem, thus 


h{t) ~ Z(o) + EC 


■Z'ih) 


( 5 ) 


If it is possible to locate the roots of the equation Z(X ) » 0, 
the expansion theorem furnishes a practical and working 
solution. In many important circuit problems, however, 
the location of the roots of Z{X) - 0, is quite impossible, 
and consequently the expansion theorem as it stands is 
only of limited practical utility! The way Heaviside over : 
came this limitation will now be shown. 

Putting t - 0, in the expansion theorem, we find that 


n 

h{a) ~W) + Ekmz) 

which represents the initial conditions in the network, that 
is, the distribution of current at the instant the switch is 
dosed. For all electric networks met with in practice 
h(o) is zero, and we may write 


1 

W) 


« 
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Substituting this value in the expansion theorem, we obtain 


h{t) = ~ Ei K z\?.' K y {l ~ eh < ] 


Expanding the exponential, h(t) becomes 
, . t yi 1 t 1 y* t* yi Ai 

() = IiZfZ'(V) + 2 iZ( r(i Jt ) t 5iLz'(y + 

1 


+ r yft: 


n\LaZ\l K y 


It is apparent that this series for h(t) may be obtained from 
the partial fraction expansion of the impedance function 
1 /Z(p) t derived in Chapter 2 : for we may write this 
expansibn in the form 




Expanding the second factor of the summation by the 
binomial theorem and arranging in inverse powers of p, 
we get 


^"pEz-^pE 


A* I v A 2 g 

Z\X K ) p^Z\l K f ■ 
i y , A n “ 1 

+ p*ZiZ'(l K ) 


+ . . . 


By comparing this expansion with the expansion in powers 
of t , we see at once that the coefficient p~ n is identical 
with the coefficient t n /n\. Consequently, it can be seen 
from these equations that if we expand the function 
1 /Z{p) in the form 


CO 



0, 1, 2, 3 . . .) 
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then the indicial admittance function is given by 


CO 



Z = 0 


which is obtained from the operational ‘expansion by 
replacing p~ n by t n /n\. 

It will be observed that in making this extension 
Heaviside discovered a powerful method of evaluating 
h(t) which does not depend upon the location of the roots 
of any equation and consequently may be used for those 
cases to which the expansion theorem cannot be applied 
readily. 


POWER SERIES SOLUTIONS 

The power series method outlined above was much used 
by Heaviside and depended upon the expansion of the 
function i /Z(p) in inverse powers of p and then replacing 
p-n by t nj n i This power series method has the great 
advantage that it can be readily carried out to obtain a 
result and involves only elementary mathematics. 

As a simple example we will now apply this method to 
a circuit consisting of an inductance L in series with a 
resistance R. The differential equation is, in the above 
notation 

L~ + Rh-[ i] 

The operational equation is obtained by replacing d/dt by p, 
and we have 

h 1 
R + pL 

This equation must now be expanded in inverse powers 
of p. Thus 


where a - R/L. 
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Writing this equation in the form 

-w 

and expanding in inverse powers of p, we obtain 

‘Wr©’*-} 

and replacing p~ n by t n /n\, we get 

which can be recognized as 

m-Ui-e-M) 


the well-known indicial admittance function for this circuit. 

If we had written the operator as 5.——^—, and then 

iv 1 + a/p 

replaced p~ l by // 1!, we should have obtained an incorrect 
result. In general, an operational equation is converted 
into an explicit power series solution only by the proper 
choice of expansion of the reciprocal of the impedance 
function. 

The process of expanding 1 /Z(p\ in such a form as to 
permit of its being converted into an explicit solution is 
what Heaviside calls ‘ algebrizing * the operator. The 
general process of expansion applicable to all cases where 
a power series solution of an electric circuit problem exists, 
may be stated formally as follows : 

Put _ 

then expanding tp(x) as a Taylor’s series, we get 
vi x ) - Vi°) + *v'(°) + j|V'"(°) + ■ • • + ^j¥’”(°) + - • • 


where 
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If now we put 

we obtain what is called an ‘asymptotic’ expansion, 
namely 

1 j 4. d. 1 4. d* 4. . d? + 

Z(f) 0 f P' p« "• 


It will be recalled that an asymptotic expansion is defined 
as the expansion of a function in inverse powers of the 
argument. The theory of asymptotic series occupies an 
important place in modern analysis. With the question of 
divergence, convergence and summability to which pure 
analysis is largely restricted we are not concerned, but only 
with Heaviside’s methods of asymptotic expansion. If the 
above asymptotic series is regarded as an expansion in the 
variable p, instead of as a purely symbolic expansion, it 
will be found to have only a limited range of convergence. 
This fact may be neglected, however, because the series 
we are really concerned with is obtained by replacing 
p~ n by t n /n\ f namely, the series 


„ A t t A t t 2 
h(t) m A 0+-J+-JT + 



which is called the ‘associated function’ of the former. 
It can be shown that if the function i/Z(p) satisfies the 
conditions for an asymptotic expansion, then the resulting 
power series solution obtained as above is always con- 
vergent.* * 

' From the foregoing example of the rise of current in 
a circuit of inductance L and resistance R, it can be seen 
that the derivation of the explicit power series is a simple 
prodess. The process of * algebrizing ’ the operators which 
arise in electric circuit theory usually involves only a straight- 
forward application of the binomial theorem. As an 
example of ‘ algebrizing ’, suppose it is required to- obtain 
an expression for the charge 0, in a series circuit f& 
resistance R, inductance L and capacity C, in response to 
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the application of a one-volt battery. The differential 
equation for the charge Q is 

L d *Q +R d Q + 2 

L di* + R -di + c ll] 

Thus replacing d/dt by the symbol />, we may write 
q(p*L+pR+^)-i, t > 0 


from which we obtain the operational equation 


Q~p*L + pR+ i,C 
which may be written as 


2 p py 


-p'l[ 1 + p + p\ 

where <j>\ = R/L, <j>z - i/LC. 

Expanding Q by the binomial theorem, we obtain 

Arranging this expression in inverse powers of p , we obtain 

0-L/-5. -f5-f2 + f2+ . . A 

Q .L\P 2 P 3 P 4 P s J 


where the a’s are readily expressed in terms of <f> x and $ 8 . 
The explicit power series solution may now be obtained 
by replacing p~ n by t n /nl This example shows that the 
Heaviside process of ‘ algebrizing ’ is formally straight- 
forward : more difficult examples are discussed in Chapter 4. 
It will be observed that in ‘algebrizing’ the original 
differential equations no consideration was given to the 
terminal conditions ; it is one of the advantages of the 
Heaviside method that the terminal conditions in linear 
networks are automatically accounted for. It ia found 
that a power series solution does not exist for all types of 
problem relating to transmission lines, although a closely 
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related series may often be derived. When, however, 
a power series solution does exist it is easier to derive than 
that of any other form of solution. A critical review of the 
scope and value of Heaviside’s power series method in 
electric circuit theory is given in Chapter 4. 

THE EXPANSION THEOREM FOR ALTERNATING VOLTAGES 
The Heaviside expansion theorem is not limited to the 
determination of the building-up of the response of a net- 
work to a direct voltage, but can be easily extended to 
include the case of the response to alternating voltage. 
Suppose that instead of applying to the network the unit 
discontinuous voltage [1] we apply the unit discontinuous 
cosine, cos eof[i]. This discontinuous cosine is illustrated 
in Fig. 2, and for purposes of analysis may be regarded 
as the real part of ^“*[1], where a)/ 271 represents the 
frequency of the applied voltage. It is shown in Chapter 8 
that if the applied* voltage is of the form £^[1], instead 
of the unit function [i], Heaviside’s expansion theorem 
takes the more general form 

nt\. Ee, 7L + fV* 
u z(jo>) +11 2j(h-j4z%) 

K *= 1 

If the applied voltage has some phase angle y), i.e. 
E cos (o)t + y>) t then the applied voltage must be written 
as the real part of the function Ee> f e imt [i\. Thus the 
expansion theorem now becomes 

... EelnW . eh 1 

■ K) “ Z{p>) +Be Z(X K -jw)Z'(X K ) 

The first term on the right-hand side of this equation 
represents the steady-state current, whilst the summation 
terms give the transient oscillations. In applying this 
theorem to the solution of any problem, it should be noted 
that only the real parts must be taken, for, in general 

' Z(jw)« R+jX 
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and taking the real part of i(t), we get 

Y v _ & cos + y> - (j>) Jy? cos (y> + 0) . e V 

VWTx* faT&Tttyi ) 

X 

where 6 - tan' 1 -^ 

K 

and 0- tan _1 ~ 

The application of this expression to a simple problem 
should make the matter clear. Consider a condenser C 
in series with a resistance R to which we apply an alternating 
voltage E sin cot, the switch being closed at the zero of 
voltage. Let it be required to determine the current 
response i(t). It can be seen at once that 



and 


Thus 

Xl ‘~CR 

Also 

• z W~Tc 


Z'lXJ'.-CR 1 


Substituting in the equation for i(t), we obtain 

E cos ^ j £e"^cos^- + ^ 

Vi?'+ i/® 2 c‘ cW<#*+i/cyF 

which reduces to 

i(t) , F.[ sin (cot - <j>) e~CR 1 

WR* + i/o)*C r a)C{R*+ i/(o 2 C 2 )) 

the well-known transient response of the condenser circuit 
to a sinusoidal voltage ^ 
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It will now be shown that we get exactly the same 
expression by using the original expansion theorem 
expressed by equation (^) and then using the h(t) function 
in the superposition theorem expressed by equation (3). 
From equation (5) it is a simple matter to show that the 
indicial admittance function h(t) is given by 

1 JL 

W-g-CR 

Now E(t) - E sin cot and E( 0) - 0, and consequently 
equation (3) may be used in the form 

t 

0 

where E'{p) - coE cos cop - ^—(e ia>p + e ®*) 

Thus we have 

t ^ 

i(t) - + e~i m p) dp 

0 

Performing the indicated integration, we find that i(t) 
reduces to the well-known expression for the alternating 
current transient as before. 

Limiting our consideration to sinusoidal voltages, it may 
be remarked that since we can obtain the formal solution 
of an alternating current transient problem directly from 
the extended form of Heaviside’s expansion theorem, why 
bother about the definite integration involved in the 
application of the superposition theorem to the indicial 
admittance function h(t) ? If the roots of the equation 
Z(A) « 0 are easily located and are finite in number and 
different from each other, then the evaluation by the 
extended form of expansion theorem is preferable. This 
expansion theorem solution, however, requires the response 
of an electrical network to be expressed as 1 sura of char* 
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acteristic oscillations or pulses, a form of solution which 
in many practical cases does not admit of easy interpretation 
or calculation, and in many trarysmission line problems 
does not exist. The superposition theorem can always be 
applied to any class of network problem provided the 
indicial admittance function h(t) has been evaluated. This 
theorem can deal with any type of applied voltage and can 
always be evaluated. If in any particular case the super- 
position integral is difficult to evaluate formally, it is always 
possible to evaluate it by graphical or mechanical processes. 
On these facts rests the practical as distinguished from the 
purely theoretical value of the superposition formula. 


SUBSIDENCE OF THE CURRENT IN A NETWORK 

The original Heaviside expansion theorem is not limited 
to the determination of the building-up of the response of 
a network to a direct voltage, but also gives correctly the 
subsidence of this response when the voltage is removed 
suddenly.; Suppose that the voltage has been applied for 
a sufficiently long period to ensure the wiping out of the 
transient oscillations, then the current in the circuit is 
given by the first term of the expansion theorem, namely 

i.JL 

1 2(o) 

Suppose that instead of removing or short circuiting the 
applied voltage, we let it remain and introduce to the 
network an exactly similar voltage but of opposite sign. 
It can be seen that the effect will be the same as if the 
original voltage was suddenly removed. This new voltage 
will produce a current given by the expression 

e K y 
‘^(0) 

Owing to the linear character of the network, the resulting 
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current will be the sum of the two currents I x and / a . 
Hence we have 



which is a formula expressing the subsidence of the current 
in the network when the voltage is suddenly removed. 



CHAPTER 4 

LADDER NETWORKS 


ARTIFICIAL LINES 

The power series method described in Chapter 3 was 
much used by Heaviside. It has the great advantage that 
it can be readily carried out to obtain a result, and the 
corresponding disadvantage that this result (unless the 
series can*be recognized and summed) is often in a form 
that converges slowly for large values of the time t. 

It is an easy matter for an engineer dealing with simple 
networks to formulate the indicial admittance as a power 
series by expanding asymptotically the reciprocal of the 
impedance function. If the network contains transformers 
or bridge sections, these should first be replaced by their 
equivalent circuits. In this chapter we propose to show 
how Heaviside dealt with such complicated structures as 
artificial lines. Artificial lines are of -practical importance 
in their use for balancing telephone cables spid duplexing 
telegraph cables and other purposes. They also find 
application in the study of electrical phenomena connected 
with transmission problems on long lines. 

The artificial line is a ladder structure consisting pf 
lumped impedance Z x *in series with the line and lumped 
impedances Z a in shunt across the line. We regard a ladder 
network as a periodic structure composed of a series of 
T sections connected in tandem as shown in Fig. 3. 
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We will assume that the voltage is applied at the middle 
of the initial or zero’th section as shown in Fig, 3. It will 
also be observed that jhe network is terminated by its 
characteristic impedance Z 0 so that the reflected wave is 
suppressed. 

In order to calculate the transient response of a ladder 
network of the type shown in Fig. 3 it is first necessary 
to derive an operational equation which formulates the 
propagation phenomena. Begin by writing down the 
voltage drop in the nth section. From Kirchhoff’s law 
we get for the nth section 

0 “ (%i + 2 ^K ” Zzih-i + hi+l) 
where i n - mesh current in the nth section 

* n -i“ » » 11 » (n - i)th section 

*n+i- » 11 ,, ,1 (n+i)th „ 

Since the network is closed by its characteristic impedance 
Z 0 , the reflected wave is suppressed, and consequently we 
may assume that a solution for i n exists in the form 
i n -ae~ ny 

where a andy are independent of n, the number of sections. 
Substituting in the Kirchhoff difference equation, we 
obtain 

. 1 ' 0- (Zj + zZJae-w - £,«-*»(«» + «-») 

From this equation it. js an easy matter to show that 

g il/ « 1 + 

where 

4^2 

In artificial line theory we are mostly concerned with the 
entering current * 0 and the current in the nth section. If 
F 0 is the voltage applied at the mid-series position of the 
zero’th section we may write 

* 0 * a m Eq/Zq 
and consequently i n - E^e “ nv /Z 0 

An expression for Z 0 may be readily obtained by considering 
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the first T section terminated by Z 0 ; we have, looking into 
the terminated T section, 


zt-\z x + 


Zz(^Zi + Z 0 ) 


( 6 ) 


Z t+ iZ 1 + Z ( 

from which we obtain 

z 0 * Vz^t + <f>) 

Consequently the expression for i n may be written as 

. = E Q {V7+4>-Vwn 
Vzm i + f) 

which is a convenient form for the expansion of an opera- 
tional equation and the investigation of transient pheno- 
mena in artificial lines. 

ARTIFICIAL TELEGRAPH CABLE 
We will now apply equation (6) to the determination of 
the indicia! admittances of an artificial telegraph cable of 
negligible inductance and leakance. Such an artificial cable 
is shown in Fig. 4, where 
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Fig. 4. — Artificial telegraph cable. 

Z x - R, the series resistance elements 


and 


Z z - the shunt capacity elements. 


In order to obtain an equation for the input current i 0 , 
we put n - 0 in equation (6), thus 

i i 0 " E 0 /V Z x Z t (i + (f>) 

Substituting for Z x and Z % and putting E 0 equal to [1], we 
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see that the indicial admittance function h Q (t) can be 
obtained from the operational equation, 




R\ 


For brevity puty - 4 /CR and then expand the denominator 
in a series of inverse powers of />, thus 


8 )" 


*/ « y + i 2 s^! + \ 

" R\ 2 'p 2*4 p* 2'4'6 ’p* * ' j 

Replacing p ~ n by t n /n\ in the above expansion, we obtain 


L2i s(y*\ , 
■ (30 2 1 


where y 1 4 /CR. 

For large values of t this series does not converge rapidly 
enough for easy computation and consequently an effort 
must be made to recognize and sum the series. At this 
point we will consider classes of functions which are defined 
by series expansions, for example 


sin x - x — r+ -7 - 


and 


a : 2 a : 3 

• I - A?+ -7 ; + 

2! 3! 


are series definitions of well-known functions. Likewise 
the Bessel functions J«(x) and I v (x) are defined, when n is 
zero or a positive integer, by the absolutely convergent 
series: 

and 

In{x)--^r-rAl + 


*» j 


2 ( 2 « + 2 ) 2 . 4 ( 2 »+ 2 )( 2 » 
* 2 X 4 


+ 4 )"‘"} 


2(2« + 2) 2 . 4(2 n + ip n + 4) 
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For large values of x, the J n (x) function is oscillatory, 
and ultimately behaves as — cos ^ The 

J n (x) function, however, increases indefinitely with x and 
behaves ultimately as tF/Vzjix. 

Now Bessel functions have been extensively calculated 
and fairly complete tables are available. Consequently, 
expressions containing Bessel functions can often be 
written down immediately by consulting tables. Thus if 
we can manage to express a slowly convergent power series 
solution in terms of Bessel functions, then the arithmetical 
work is greatly facilitated. The series expansion given for 
h 0 , however, is so slowly convergent for very large values 
of t that its use for calculation purposes is not justified and 
consequently we must try to express this series in terms 
of tabulated functions. 

If we put n m 0 in the series expansion of I n (x) } and then 
multiply I Q (x) by e ~ x term by term, we obtain on combining 
terms of the same powers of x , the expansion 




LU X 

( 30 2 


which can be made identical with the series given for h 0 
by putting x - yt/ 2. Thus we may write for the power 
series expansion of A 0 , the compact formula, 


2 


which can be rapidly evaluated for all values of t from the 
tabulated values of the function e~ x I 0 (x) given in books 
of mathematical tables. • 

We will now derive the power series solution for the 
indidal admittance of the nth section. We begin by 
writing equation ( 6 ) in the form 


MO- 


(Vi+^+ v^)2 s " 


from which we will develop the operational formula. 


4 
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Now (Vi + <£ + V<f)~ in - <f>~ n (i+ V i + <j > ~ i y 2n 
and VZA(i+ = V^ZjZjV i + (f>~* 

Since V <f>Z x Z % » .ft/2 

We may write 

where the factors and ft 2 (/>) are defined as 


h ' {f) V [ + 4/pCR 


and 


FJj>)-(T + Vr+4/pCR)-*’' 


In order to derive the required power series expansion 
each factor must be expanded by the binomial theorem and 
arranged in inverse powers of p. After this the two 
resulting series must be multiplied together and rearranged 
in powers of l/p. Applying this process, we get the series 
, 2 (7 2 \ w (zn + z)/ 2 \ n + l 

n ' m zW\[pCR) ~W[pCRj 

(zn + 3)(2n + 4)/ 2 W 2 

+ ^2!) [jCR) 

Replacing p~ n by t n /n\ } we get 

2 fl/2*\ n (2tl+2) /2f\ n41 

w " 2^{id\Cft ) ~ 2(i!)(» + I)! \^Cft ) 

(2»4-3) (2» + 4)/ 2 t\ nVi 

z\z\)(n + 2)! \CrJ 

For large values of / and n this series is very difficult to 
compute or interpret, but by a similar process tor the one 
outlined for the input indicial admittance the above 
power series can be recognized as equivalent to the series 
expansion of the function e~ X l n (x). Hence we may write 
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where I n (zt/CR) is the modified Bessel function of the 
first kind of order n and argument (zt/CR). By putting 
n o, the above equation reduces to h 0 (t) as it should 
do. The input current i 0 (t) and the output current i n (t) 
in response to the direct voltage £ 0 shown in Fig. 4, may 
be obtained from h 0 (t) and h n (t) by multiplying these 
expressions by the constant E 0 . If, however, E 0 is an 
arbitrary function of the time, these current responses can 
be obtained by substituting h n (t) in the superposition 
theorem expressed by equation (3) and performing the 
indicated integrations. 

ARTIFICIAL COIL-LOADED TELEPHONE LINE 
Suppose now that the resistances which constitute the 
series elements Z x of the ladder network are replaced by 
.inductances as shown in Fig. 5. 



This type of artificial line is closely related, in its 
properties and performance, to the periodically loaded 
telephone line, and its solution is, to a first approximation, 
a working solution for the coil-loaded line invented by 
Heaviside. 

In this case we have for the operational impedances 
Z x *= R + pL 

and Z a “ 1 /pC 

The power series solution. for the current response i n in 
the nth section may be developed from equation (6) as 
before. This equation may be written in the form 
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Hence we get for the denominator 

i /CW i h 

Now 

, / / /coA 2 i \ _2w 

(j + + + *J 1 + \j) {i+R/pL)) 

Hence the equation for i n may now be written as 



which is a convenient form for asymptotic expansion in 
inverse powers of p. In order to derive the required 
power series solution each of the last three factors must be 
expanded by the binomial theorem and arranged in inverse 
powers of p . The next step is to multiply the three result- 
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ing series together and rearrange in inverse powers of p. 
This series must then be multiplied throughout by the 
factor (<o c /P ) 2n+1 an d then the explicit power series solution 
may be obtained by replacing p~ n by t n /nl This expansion 
becomes very complicated if a number of terms is required, 
and it will be found that owing to its slow convergence 
the series is practically useless to compute from except for 
small values of the time t associated with the first transient 
surge. Consequently, it becomes necessary to make an 
effort to recognize and sum the series, which in this case 
is a most difficult matter. It will be found that t n can be 
written as 

X 

- B, Jc e 

0 

where fjy/l 

X - (O c t 

and J and I are Bessel functions. The evaluation of this 
expression depends on the numerical computation of Bessel 
functions and the integration of definite integrals. The 
latter operation presents no difficulties and the integrations 
are easily and accurately effected, once the integrands are 
plotted as time functions, by means of a suitable integraph 
(or even on the drawing-board). Fortunately, also, the 
Bessel functions have been the subject of exhaustive study 
and tabulation. The formal solution of the integral would 
consist of expanding the Bessel functions as power series, 
integrating term by term and rearranging. This would 
yield the formal Heaviside power series solution we have 
discussed above. It should be noted that the expression 
for % is a special case of a more general formula deduced 
by Dr. J. R. Carson in his researches on speech-excited 
transients in transmission networks. For a complete 
discussion of the artificial coil-loaded line the student is 
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referred to Carson’s paper on ‘ Transient Oscillations 
Transactions A.LE.E . * 1919. 

The Heaviside operational scheme we have outlined in 
this chapter . has several distinct advantage's from an 
engineer’s point of view. It can be readily carried out to 
obtain an explicit solution. It involves only elementary 
mathematics. In many cases the results obtained would 
be very difficult (or even impossible) to obtain by iny other 
method. The scheme, however, has certain disadvantages. 
The results obtained are sometimes in a form that is not 
convenient for numerical study owing to the slow con- 
vergence of the series ; if, however, the series can be 
summed in terms of tabulated functions, this difficulty can 
be avoided. 

It must be noted that a power series solution does not 
exist for all classes of electric circuit problem ; for example, 
in a certain class of network problem relating to transmission 
lines, a power series solution does not exist. In this class, 
however, a closely related series in fractional powers of the 
time t can usually be derived. A discussion of the applica- 
tion of this related series to network problems will be 
found in Chapter 7. 



CHAPTER ^ 

HEAVISIDE’S LAST THEOREM 

IMPULSE FUNCTIONS 

In his researches Heaviside did not bother to construct 
rigorous mathematical proofs of his theorems ; to him a 
‘ rigorous proof * was only an attempt to meet the whims 
and fancies of certain pure mathematicians. Heaviside 
‘ proved ’ nothing ; he simply made sure he was right. 
To make sure of being right he appealed to experiment 
and measurement for a judgement : these practical tests 
he applied to his theorems in many ways and under many 
different conditions. When dealing with his asymptotic 
series, for example, he depended upon a practical grasp of 
the physical side of the problem he was analysing to guide 
him away from mathematical pitfalls. In many cases he 
used physical arguments to establish mathematical results. 
Thus conflict with pure mathematicians was inevitable. 
This conflict was unfortunate because it antagonized 
him, and consequently he made no attempt to pre- 
sent his operational calculus as a rigorous mathematical 
structure. 

In this chapter attention is drawn to an important and 
far-reaching Heaviside theorem which so far appears to 
have escaped the notice of engineers. This theorem has 
been reconstructed from the scattered papers of Heaviside 
and is probably the last theorem he ever deduced. For 
this reason it has been called ‘ Heaviside’s Last Theorem \ 
From the basis of this theorem it is possible to deduce 
rigorously all' the various Heaviside operational processes 
used in electric circuit theory. 

At this stage it would be well to consider the nature of 
the* operator p ( or d/dt) which is assumed to act on real 
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variables. The inverse operator p~ l represents the opera- 
tion Jdf. It follows, therefore, that 

i 

o 

and consequently the operator p undoes the operation ^ -1 
if p .acts after p~ l . When p and p~ l both occur in an 
operator the integrations must be carried out before the 
differentiations. 

t 

Otherwise />->/>//(«) = Jtf'(t) dt - H(t) - H(o) 

0 

which shows that p and p~ l are commutative only if .the 
function operated on vanishes with t, i.e. H( o) - o. It 
will be noted that Heaviside’s unit function [i] preserves 
this commutative property in the operators which arise 
in electric circuit theory, because 

[i] * i for t > ol 
* o for / < oj 

Thus [i] is zero everywhere in the range, - oo < x < -f oo, 

except at the point x - o ; at this point the unit function [i] 
suddenly jumps from zero to unit value as shown in Fig. 2. 
The derivative of the function is therefore infinite at the 
point t * o. From the above 

/>[i] « o for t >-o 

- oo for / - o 

But rWiD-rW-ff" 

Thus p[i] is a function which is everywhere zero, expept * 
at the origin where it is infinite ; the integral of the function, 
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however, is unity. This function, which Heaviside called 
a ‘ unit impulse is illustrated in Fig. 6. 

To visualize the unit impulse function consider a long 
thin strip of length / and width dt extending along the 
axis of y as shown in Fig. 6. Now let / — ► oo as dt — >• o 



Fig. 6. — The unit impulse. 

in such a way that lim (/ dt) « i. From Fig. 6 it can be 
seen that the derivative p[i] is composed of three sections : 

(1) The f-axis from / * - co to / = o. 

(2) The limiting form of the strip shown, where / — ► 00 

and dt —> 0 in such a way that the area (/ dt) — ► 1. 
This strip is perpendicular to the f-axis from the 
point t « 0. 

(3) The f-axis from / - 0 to t » + 00. 

This mathematical idea may be illustrated by considering 
the charge that passes to the plates of a pure condenser 
when connected to an e.m.f. Thus if q is the charge, 
i the current and e the voltage, then 

e 0 
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If C and e have unit values, it follows that 

Thus the current p [i] represents a unit impulse, i.e. it is 
infinite at time t * o and zero thereafter. This means that 
at time / - o, when the switch is closed, a unit charge is 
instantly placed on the plates of the pure condenser. 

TRANSFER OPERATORS 

The symbolic form of Taylor’s theorem was well known 
in Heaviside’s time. It is easy to show that f (x± vt) is 

re- 
produced by the operator^ e dx operating on f(x). 

Thus e ±vt dx f(x) « f(x ± vt) 

If t is the time and v a velocity, then the expression 
fix - vt) represents a curve or wave moving along the axis 
of x with a velocity v in the positive x direction. In general, 
the effect of the operator e ~ ph on fit) is to translate it 
bodily through a distance h to the right. 

Exponential operators of this type are known as transfer 
operators and allow functions to be moved about in time 
or position. Heaviside made great use of these operators 
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in the practical applications of his Duplex electromagnetic 
theory. He used them to remove damping factors; to 
transform space waves into travelling waves ; to make 
travelling waves stand still so that he could study their 
distortion with lapse of time. 

Now consider the function e~ p \i ] ; the effect of the 
operator e~ ph on [i] is to move the unit function bodily 
through a distance h to the right as illustrated in Fig. 7. 
This transferred unit function can be written [i] ft ; then 

Likewise the effect of the operator e l>h on p[i] is to move 
the unit impulse from the origin to a point on the /-axis 
where / «= h. In other words, the effect of the operator is 
to transfer p[i] with all its properties unchanged to another 
origin of co-ordinates at the point ( h , 0) as shown in Fig. 7. 


SPOTTING FUNCTIONS 


From Fig. 7 it can be seen that the transfer operator 
e ph applied to the unit impulse p[i] generates a function 
composed of three sections : 

(1) The /-axis from / = - 00 to / - +h. 

(2) The limiting form of the unit area strip extending 

perpendicular to the /-axis from the point / - + h. 

(3) The /-axis from / * + h to / - + 00. 

Heaviside called this function a ‘ Spotting function ’ 
because it picks out or spots a certain value of a function. 
In order to illustrate this idea, begin by considering two 
real variables, such as electric and magnetic forces, which 
range from 0 to 00. For simplicity denote these variables 
by x and /. Let a denote d/dx and p denote d/dt. Then 
if [1] represents the unit function, we have 


<x[i] - unit impulse at the point x « 0 
p[i]~ » „ „ „ „ *-0 

e'M i]~ „ „ „ „ „ x-t 

» >• •• » 


»» 1» 
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It is apparent that i] is equal to e~ px p[ i], because 
a unit impulse placed at the point x » t on the x axis is 
the same thing as a unit impulse placed at a point x 
on the t axis. Thus e~ d a[i] is a function of x concentrated 
at the point x - t ; while e~ px p[i] is a function of t con- 
centrated at the point t - x. This fact makes the two 
spotting functions quantitatively equal ; Heaviside demon- 
strated this by showing that both forms of spotting function 
can be obtained from the same operational generator, 
op/(o + p). By expanding this operator in inverse powers 
of p and then replacing p~ n by t n /n\, the spotting function 
e - j 8 produced. Likewise, by expanding the operator 
in inverse powers of a and then replacing a~ n by x n /n\, 
the spotting function e~ px p[i] is produced. Thus both 
spotting functions are generated by the same operator and 
are thus quantitatively equal 

= «' "M 1 ! - e ~ rx P [i] 

From the above it can be seen that a unit impulse is a 
function which is everywhere zero, except at a single point 
where it « infinite. The integral of this function, however, 
is unity. 

AN INTEGRAL THEOREM 

Consider two real variables x and t ranging from o to 
+ oo. Then /(#)e” rt </[i] is zero over all the range 
o < x< + oo, except at the point x-t. At this point 

00 

f(x)e ~ ^[i] is indeterminate. The integral j/(#)e ~ d a[i]dx 

, - 00 

can, however, be evaluated. From Fig. 8 it can be seen 
. that there is no contribution from the range o < x < *, 
or from the range t < x < + oo, since the integrand is 
zero over this range. To obtain the contribution at . the 
point t, consider the function as a thin strip 

of length / and width 6x, The area of the strip under the 
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integrand /(*)*“ *<t[ i] at the point x**t, is therefore of 
width 6x and length l f(t) as illustrated in Fig. 8. In the 



Fig. 8. — The function y —f(x)e~ot a[i], 

limit / — > oo as (5 a; — > o, in such a way that lim (/ <5#) — ► i , 
and lim {/(/)/ d*} — > f(t ) ; thus we may write 

00 

O' 

which is apparent when we observe that the function to 
be integrated is an impulse of size f(x) at the point x-t, 
at which point f(x) has the value f(t). This leads us to 
Heaviside’s last theorem which may be stated formally in 
the following manner. If x and t are two real variables 
ranging from o to oo, and a denotes d/dx and p denotes 
d/dt, then 



It will now be shown that Heaviside’s last theorem is a 
fundamental one, i.e. the ideas involved in it are related 
in a natural way to a large complex of other mathematical 


54 Heaviside’s electric circuit theory 

ideas. Thus the theorem leads to'many others ; not only 
electric circuit theorems, but also to theorems involving 
Fourier and Bessel integrals, elliptic functions, &c. In 
this book, however, we are only interested in its applications 
to electric circuit theory. 

Since Heaviside’s time the two most original contribu- 
tions to electric circuit theory were made by Professor 
I’ A. 1 Bromwich of England and Dr. J. R. Carson of 
America. Bromwich’s contribution will be found in his 
paper, ‘Normal Co-ordinates in Dynamical Systems’, 
read to the London Mathematical Society in 1916. In 
this paper Bromwich deduces an integral theorem from 
which it is possible to formulate an explicit relationship 
between the indicial admittance function and the impedance 
function. Bromwich’s discussion of Heaviside’s methods 
by means of contour integrals in the complex plane is one 
of great beauty and power. It justifies and extends many 
of Heaviside’s processes. Nevertheless, it would be 
difficult to imagine two methods more unlike than those 
of Bromwich and Heaviside. In this connexion it is 
interesting to note that Bromwich has repeatedly declared 
that the direct operational method of Heaviside is far the 
best for dealing with the class of problems concerned. 
Carson’s method, on the other hand, is very similar in form 
to that of Heaviside. It approaches the subject in an 
original manner, and, by basing the analysis on what is 
known as an integral equation it puts to new use a largo- 
number of operational results in such a way as to greatly 
extend their range of application. Carson’s contribution 
to electric circuit theory will be found in his book, Electric 
Circuit Theory and the Operational Calculus , published in 
1926. In this book Carson deduces an infinite integral 
theorem which expresses an implicit relationship between 
the indicial admittance function and the impedance func- 
tion. From this theorem Carson establishes the operational 
formulae given by Heaviside and presents the operational 
calcuhis as a rigorous structure. In his book, however, 
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Carson has failed to notice that his fundamental integral 
equation is only a special case of Heaviside’s last theorem 
and that Bromwich’s contour integral is the solution of 
his equation. 

Heaviside gave the Carson theorem on page 236 of the 
third volume of his Electromagnetic Theory , published in 
1912 ; his notation and discussion at this point, however, 
was so different from his usual scheme that the theorem is 
not apparent until it has been transformed into a more 
familiar notation. Heaviside did not use the theorem 
extensively for electric circuit analysis ; unlike Carson, he 
did not use it to establish or transform operational formulae, 
but used it to evaluate infinite integrals. This may be due 
to the fact that he deduced the theorem rather late in life, 
or he may have preferred other methods. In this con- 
nexion it must be noted that there are many perfectly valid 
methods of establishing operational formulae without using 
the integral theorem, for example, the expansion theorem 
and power series methods discussed in Chapters 2 and 3. 
Consequently the use Heaviside made of the theorem, 
namely, to evaluate the infinite integrals which arise in 
electromagnetic theory and which cannot be evaluated 
readily by any other means, may become the more valuable 
process as the operational calculus develops. In any case 
the theorem is of use either way. Since our subject is the 
analysis of electric circuits, and not the solution of infinite 
integrals, we return to a discussion of Carson’s integral 
theorem. 


carson’s integral equation 

Applying equation (7) to electric circuit theory, we 
obtain ' 

CO 

I(J)~\ KX)e ~ p ’ fdx 

0 

by replacing h(t) by its generating impedance function 



56 Heaviside’s electric circuit theory 

i/Z(p). Since the integration is with respect to x, the 
operator p that occurs can be treated exactly as though it 
were a parameter; and in fact it can be replaced by a 
parameter A. Thus we obtain on replacing x by t 

00 

h(t)e-»dt (8) 

0 

which is Carson’s integral equation. It can be seen that 
this equation expresses an implicit relationship between 
the impedance function and the indicial admittance func- 
tion h(t). By virtue of this relationship Carson puts to 
new use a large number of established operational solutions 
in such a way as greatly to extend their range of application. 
Moreover, Carson introduced Borel’s theorem into the 
analysis and thereby gave the method much additional 
power. 

The formal solution of Carson’s integral equation is 
given by the Bromwich integral, where 

* m C+ju 0 



C -fco 


which expresses an explicit relationship between the indicial 
admittance and the impedance function. The Bromwich 
solution is discussed in detail in Chapter 8. For the present 
it can be noted that the above equations are valid for all 
values of A, for which the real part is greater than some 
finite constant C. This constant C must be at least large 
enough to make the infinite integrals converge. For the 
majority of passive networks this constant may be taken 
as 0 ; if the network employs some form of feedback, 
however, the above equations are valid only when C is 
greater than some finite value. 

The foregoing analysis naturally suggests that the 
methods of Carson and Bromwich are closely related -to 
Heaviside’s last theorem. It also follows that a Heaviside 
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operational equation is a shorthand way of writing down 
infinite integral equations or contour integrals. If this fact 
is kept in mind the operator p ^ill lose its mysterious 
character: thus if we try to evaluate the operational 
equation Z(p) = V/>[i], it does not mean that we are 
seeking to find the square root of the differential coefficient 
of a discontinuous function, which is meaningless, but 
simply that we have written down in shorthand the Carson 
integral equation, 

00 

-L -*<//, 

0 


or the Bromwich contour integral 

C+joo 


m 


i 

27 lj 


C -j'oo 


xVx 


= dX 


which shows that the operator p has lost its original 
significance and is now the complex argument of functions 
which obey all the laws of mathematical analysis. 


5 



tHAPTER 6 


THE ESTABLISHMENT OF OPERATIONAL 
FORMULAE 

SURVEY OF FUNDAMENTAL FORMULAE 

In Chapter i the analysis of a circuit problem was defined 
as the determination of the variation of a physical magni- 
tude in terms of time and some dimension ; more precisely, 
the determination of the oscillations of a linearly connected 
system described by a set’ of differential equations with 
constant coefficients or a partial differential equation of 
the wave motion type. It was shown that these descriptive 
differential equations may be readily written down by 
an application of KirchhofTs laws to the problem. The 
electric network considered was assumed to be in a passive 
state when at reference time t - o, a unit electromotive 
force [i], was suddenly applied. The first step in the 
operational analysis of a network problem is the purely 
formal one of replacing the differential operator d/dt by 
the symbol />, and thereby reducing the differential equa- 
tions to an algebraic form. From the rules of algebra we 
obtain the Heaviside operational equation which may be 
written formally as 

(+) 

This equation is a purely symbolic expression and its 
interpretation means a transformation to yield a function 
of the time t no longer involving the differential operator^. 
In this expression h(t) represents the indicial admittance 
function, while Z(J>) is the impedance function and is the 
means whereby the parameters and the connexions of the 
network enter the mathematical formulation. Thus Z{p) 
completely specifies the network ; if p is replaced by jeo, 
5 » 
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where co/zn is the frequency, then Z(jco) is the well-known 
steady-state impedance function of alternating current 
theory. From Chapter i it will be recalled that the 
Heaviside unit discontinuous function [i] should appear 
as the operand on both sides of equation (4) ; but since 
this unit function, is always tacitly assumed to apply to 
all the operational equations which arise in electric circuit 
theory it need not be explicitly stated. 

In Chapter 5 it was pointed out that an operational 
equation is not a true algebraic equation and has no literal 
meaning by itself ; it is simply an easy shorthand way of 
writing down a set of descriptive differential equations 
which have been obtained by applying physical laws to the 
network. It was also shown as a corollary to Heaviside’s 
last theorem that to every operational equation (4) there 
corresponds the Carson integral equation 

00 

Wf x \ eUh ^ dl • ■ • ■ ( 8 ) 

0 

which expresses an implicit relationship between h(t) and 
Z(p). By virtue of this relationship Carson has shown how 
to extend the range of operational equations and ratify 
unestablished operational results. 

The formal solution of Carson’s integral equation is 
given by the Bromwich contour integral, where 


C+joo 



C-j 00 


which expresses an explicit relationship between h(t) and 
Z(p). By virtue of this relationship Bromwich has shown 
how to evaluate Heaviside operational equations by direct 
integration and thereby determine h(t). 

Having determined h(t), the response i(t) of the network 
to any arbitrary applied electromotive force E(t) is obtained 



from an application of the superposition theorem 
t 

'(>) = - #M dr • • • (3) 

0 

These equations constitute a complete mathematical formu- 
lation of the problem of electric circuit analysis ; and their 
solution automatically takes care of the necessary terminal 
conditions. That is to say, contrary to the classic method 
of solution discussed in Chapter 1, the explicit determina- 
tion of roots of high degree algebraic equations and the 
explicit determination of integration constants is rendered 
unnecessary. 

It will now be apparent the expansion theorem and power 
series solution discussed in Chapters 2, 3 and 4 are examples 
of the fact that if an operational equation can be expanded 
in a sum of terms, thus 

111 1 

m~m + zjp) + "' + ^) 

and if the auxiliary operational equations 

**■ w K ~ utv ' m 

can be solved, then 

K { ) -*! + *> + *« + ■ ■ ■ +K 

Thus we can juggle the left-hand side of Carson’s integral 
equation (8) in any way we please provided the changed 
expression is actually equivalent to 1 /Z(X). We may 
develop in series, split into partial fractions or the like, 
all without altering the validity of Carson’s integral equa- 
tion. Thus an operator may be expanded into a sum of 
terms and to each separate term will correspond a valid 
integral of the Carson type. 

SHORT TABLE OF CARSON INTEGRAL EQUATIONS 
A large number of infinite integrals of the Carson type 
have been worked out by mathematicians. Consequently 
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the solutions of many of the operational equations which 
arise in electric circuit theory can frequently be written 
down by inspection. When this h not the case, however, 
the proper procedure is then to expand the function 
i /Z(A) by any means, in such a form that the individual 
terms are recognizable as being symbolically equivalent 
with Carson integral equations. It must be noted that 
not all the operational equations that arise in electric circuit 
theory can be recognized without transformation ; this 
idea corresponds with the fact that a table of integrals is 
not always sufficient to recognize an integral, but must be 
supplemented by general methods of transformation and 
integration. Below is a small list of some Carson integral 
equations. This list could be greatly extended. Even 
this short list, however, contains the solutions of many 
operational equations which occur in important technical 
problems. 


A 

(X + by 




e hi .dt. 


/? -2 
VX.e =A 


.dt . 


J)/t 

Vjrt 


(*) 


'’(y^J-X^-«J„(k)dt . . . ( c ) 

0 

where r 2 = A 2 + b 2 


A 

VFTi 
where / > b 


e-W'*' ~Xy-uj 0 {Vt 2 -b 2 )dt . (D) 
o 


00 

r^-^e-».e-u m dt ...(E) 

V T « 

where J n and I n are Bessel functions and b a real quantity. 
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DEVELOPMENT OF INTEGRAL SOLUTIONS BY OPERATIONAL 
METHODS 


We will now discuss in greater detail how Carson’s 
integral solutions may be built up from the operational 
methods discussed in Chapters 2 and 3. Consider the first 
solution (A ) ; putting n « 0, we have 

an 

X+b J 
0 

from which we obtain the operational solution 


_ e ~bt 

p + b 

which can be easily verified by an application of Heaviside’s 
expansion theorem. In order to obtain the operational 

solution of ~~~T , from the above equation, write --^-7 
P + b* ^ p + b 


p'p + h ‘ 
solution 


j. We have by associating p 1 with j dt, the 


1 ~ [e MJt- !(i - e w ) 

P+b J b x 
0 

This operational solution is a well-known result, but the 
method of writing the operational equation as the product 
of two operators and then evaluating should be carefully 
noted. 

Suppose we transform an operational equation into the 

form of the product of two operators, namely, 

and vrite the results of the individual operators as 
/»!(/)- i/Zi(j>) and h 2 (t) ■= i/Z 2 (p). Then we may write 
the operational product as 
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Now this equation may be regarded as the effect of suddenly 
applying an electromotive force h 2 (t) to a network the 
impedance function of which is 2^(p). Consequently the 
response of the network may be obtained from the super- 
position theorem expressed by equation ( 3 ). Thus we 
may write 

t 

zjf) ' = 7t\ _ v) d v 

0 

It will be noted that this expression involves a variable 
parameter under the sign of integration and that ip is the 
variable of integration which disappears when the integra- 
tion is performed. The expansion and evaluation of this 
type of integral is discussed in Chapter 8 where it will be 
shown that the integral expression above can be written 
in the expanded form 

t 

ZjpyW) ~ + J*i(* - 

0 

Since we have evaluated i/(p + b) and p/(p + b) } it will be 
observed that this expansion will enable us to evaluate 
operational expressions like p/(p+ b) n+1 , where n is any 
positive integer, and consequently establish Carson’s 
integral equation (A). 

Let n» 1 , then we have 

P r 1 P__ 

(p+b) 1 p + b'p + b 

where - ^(i - e~ M ) and h t (t) - e~ il . Inserting these 

results in the expansion of h 2 (t)/Z(p), and performing the 
indicated integrations, we find that 


P 

(p+by 


te~ bt 


which agrees with the Carson equation (4) for n » I. 


We 
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can continue the process by evaluating the operational 
expression p/(p + b) z , and show that we get the same result 
as would be obtained * by putting n «= 2, in the integral 
equation (A). Writing the operational expression in the 
form 

J_ _ ± 

p + b [p+ by 

we have h^t) - ^(1 -- e~ bt ) and h^t) ~ 1 e~ bi 

Inserting these results in the expansion of hJfy/ZAp), and 
performing the indicated operations, we find that 

. P -tlg-bt 

(p + by a! 

which agrees with the result obtained from equation (A) 
when n - 2. 

Continuing this process, we obtain the general operational 
solution 

(p+b)* (n-i)f 

If the operator p is replaced by the parameter A, this 
solution may be written as the infinite integral 

(A + ft)”+ 1 J n\ 

0 

which is a convenient form for transformation. Putting 
b ■ 0, in this expression, we get the well-known result 



The above shows how the product of two operators may 
be transformed by means of the superposition theorem. It 
must be noted, however, that it is possible to make these 
transformations in a much easier manner. 
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borel’s theorem 

There is a theorem due to Borel which is more con- 
venient for dealing with the prdduct of operators than 
the Superposition theorem. We begin our discussion of 
Borel’s theorem by expressing it in the terminology of 
this book. 

Suppose that the function i/Z(p) can be transformed 
into the form 

111 1 

where h x (t) - i/Z x (p), h 2 (t) = i/Z 2 (p), and p 1 denotes the 
t 

operation j dt. 

0 

Then from Borel’s theorem we may write 

t 

h(t) = j j\(# 2 (* - V’) d 'P 
0 

t 

= j*i(* - y))h 2 (y>) dip 
0 

where ip is the variable of integration. 

To illustrate the use of this theorem, let us evaluate 
again the operational expression p/(p + by. In this case 
we write 

i f _P_ P__ 

(p + by ~p'p + b (p + by 

where h x (t) ~ e~ bi and h 2 (t) = te~ u 
Substituting in Borel’s theorem, we get 
t 

h(t) - bv (t - ip)e~ b{i ~ v) dip 
0 

which reduces to t 2 e~ bl /2\, as before. It should be noted 
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that it is' much simpler to use this method than the 
superposition theorem. 

We have discussed in some detail certain operational 
processes which lead to the formulation of integral equations 
of the Carson type. These infinite integrals should be 
carefully tabulated. Thus when a new operational expres- 
sion is met it may be interpreted at once by reference to 
the table. If, however, this cannot be done directly, it 
will be necessary to pick the most likely-looking infinite 
integral and transform it into the required form. An 
illustration of the methods of transforming these infinite 
integrals will be given in the next chapter. 


Heaviside’s ‘shifting’ transformation 


In order to reduce the work involved in expanding 
operational expressions containing exponential damping or 
attenuation factors, Heaviside introduced what he called 
a ‘ shifting ’ transformation. This transformation con- 
sisted in making certain changes in the operator and 
writing the exponential factor outside; these changes 
generally resulted in simplifying the operator and thereby 
making it easier to interpret. We will describe this process 
by means of the Carson integral. It has been shown that 
if the operational equation 


exists, then 


m- 


m 


00 


i{iy 



it 


where A is a parameter. Now replace A by (A + a), we get 


0 
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which may be written in the form 

00 

a'Z(I7^j = dt 

I + /t " 

from which we obtain the operational equation 


e ■ «%(*) = 


.A 1 

/> + a’Z(/> + a) 


Replacing h(t) by its operational generator i/Z(p), and 
changing the sign of a, we obtain 


JL = e -at _J !_ — — 1 — (jo) 

Z(p) p- v. ' Z(p~v)' • • • (I ° ; 

This transformation enables us, when we have an opera- 
tional expression, one factor of the result of which we 
suspect to be exponential, to write down that part immedi- 
ately and make the corresponding changes in the operational 
equation, formulae for shifting other functions may also 
be developed. 

As a simple example illustrating the application of this 
exponential transformation, consider a circuit consisting of 
a resistance R , an inductance L and a capacity C in series, 
to which the unit electromotive force [1] is suddenly 
applied. The operational equation is 

R + pL + i/pC 

___ P 

■ifr-JWfe-S) 

where ^i = -a+//? 

A 8 =-a -jP 

and a = R/iL 

o /V~JP 

p y Z"c 4Z, 2 
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Thus h(t) may be written in the form 

wj.i Jl 

w pL ip+zy+f) 1 

We suspect the presence of the factor e ~ cd , and using the 
Heaviside shifting transformation, we write 


' pL p -a (p - a + a ) 2 + p 2 


h(t)° 


pi/ 


, pl 


It is a simple matter to evaluate the operational part o! 
this expression. Expanding in inverse powers of p and 
then replacing p ~ n by t n /n \ , we find that 



Thus we may write 


h(t)~ p.e^mpt 


This result is for the case where > a 2 . If, however, 

a 2 > we note that P now becomes a pure imaginary. 

Consequently we put P « jp i} in the above equation, and 
obtain 

*(<)-" J-jf sin j/M 

Since sin jp x t ■* j sinh p x t , we obtain 

h(t)--~e~« t sinh p x t 

These results are in agreement with those obtained by the 
expansion theorem in Chapter 2. The ease with which 
the shift transformation enables these expressions to be 
arrived at should be noted. 
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FRACTIONAL ORDER DERIVATIVES 

In Chapter 3 it was shown that the initial response of an 
electric network may be obtained by putting t « 0 in the 
indicial admittance function h(t). Consequently the func- 
tion h( 0) represents the distribution of current in the 
network at the instant the switch is closed. For all practical 
networks met with in practice h( 0) - 0, because inductances 
will carry no current and condensers will support no 
voltage until a finite time has elapsed. The same effect 
may be obtained by considering the network response to 
an infinite frequency because for this limiting frequency 
the inductances will carry no current and the condensers 
will support no voltage. Thus if Z(oo) represents the 
steady-state impedance function for an infinite frequency, 
we may write 

zk = h{0) = 0 

because for all practical networks Z( 00) - 00. 

Let us now consider some case in which Z(oo) is either 
finite or zero. For a pure condenser, Z(jco) » 1 /jouC, and 
for an infinite frequency Z( 00) « 0. Thus it is only possible 
to obtain such a result by applying a voltage to a condenser 
which has neither resistance nor inductance. Since it is 
impossible to build a condenser entirely devoid of resistance 
and inductance, it follows that the case where Z( 00) comes 
out zero is a mathematical fiction. No matter how small 
this inherent resistance and inductance may be, it will 
always be sufficient to eliminate the difficulties associated 
with infinite currents at time t =* 0. This follows because 
R +ja)L must always be added to 1 /jcoC, and consequently 
Z(oo) is now infinite, however small R and L can be made. 

Because the condition where h( 0) = 00 cannot exist 
physically, it must not be thought that this case is without 
theoretical importance. It will be recalled that in Chapter 5 
we showed how Heaviside reasoning in terms of impulses 
associated with infinite initial currents established many 
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important mathematical results. We discussed the case of 
impulsive response in a lumped system in which Z(oo) = o. 
We proceed now to the case of impulsive response in a 
distributed system with the object of introducing the reader 
to Heaviside’s fractional order derivative operators. Let it 
be required to determine the current i 0 entering a non- 
inductive cable. Now a non-inductive cable has a dis- 
tributed series resistance R and shunt capacity C per unit 
length, F 0 /Z 0 , where V 0 is the applied voltage and 
Z 0 the characteristic impedance. This impedance is given 
by the expression VR/jcoC. In order to obtain the opera- 
tional equation, put V 0 equal to the unit function [i] and 
replace jo by p in Z 0 . By making C numerically equal 
to R } we obtain the operational equation 

h{t) - Vp[i] 

The question now arises as to what meaning can be given 
to this operational expression. We have shown that an 
operational equation is a shorthand statement of an analytic 
equation. In this analytic equation the operator p has lost 
its original significance and has become the complex 
argument A of a function obeying all the laws of analysis. 
Hence from Carson’s integral theorem we may write 
00 

= ll h(t) it 
0 

The value of h(t) may be determined in a number of 
different ways. Perhaps the simplest is to refer to a table 
of Carson integrals and pick out the solution. From the 
short table given in this chapter we cannot find the exact 
form, but formula ( B ) resembles it. If we put b - o in 
formula ( B ), we have the required form, namely 


oo 



0 
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from which we obtain the result 




L 

Vnt 


This result is an important one and may be easily verified 
by direct integration of the Carson integral. By putting 
t = co 2 in this integral we obtain 

00 

VI « do 

0 

Now put x**coVI, and we obtain 
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e~ x} dx « — , it can be 

0 

seen that the Carson integral reduces to an identity. Thus 
the interpretation of Vp as 1 /Vnt is rigorously established. 

From this interpretation of Vp it can be seen that the 
entering current h(t) has an infinite value at time t « 0, 
and then dies away in accordance with the formula 1 /Vnt. 
This infinite initial current results from our ignoring the 
distributed inductance of the cable, which, no matter how 
small, keeps the initial input current finite. Now the 
charge q 0 that enters this cable may be readily obtained by 
operating on h(t) by p~\ 

We have 


P Vnt 
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i 

Associating/) 1 with | dt we obtain the explicit solution, 
°* 

o 

It can be seen from this expression that the charge that 
enters the cable is finite for all positive values of t , although 
h(t) is infinite at time t - o. Thus the operator for q 0 is 

-V p or i/Vp, which must be interpreted as 2 ^Jt. 

When Heaviside encountered operational expressions 
like Vp for the first time in his transmission researches, 
his method of interpretation was very different from the 
one outlined above. In his early researches Heaviside 
made no pretence to rigorous formulation and allowed only 
physical and vectorial ideas to guide his mathematics. 
One of his methods was to obtain the explicit solution of 
a number of problems by classical methods and then 
compare them directly with their operational equations. 
Thus he was led to attach a definite meaning to his 
operators, which significance he would apply to those 
problems in which the classical methods failed him. An 
example will illustrate this kind of inductive reasoning. 
The artificial non-inductive telegraph cable discussed in 
Chapter 4 can be made to simulate the corresponding 
smooth cable by making the number of sections n very 
large and the C and R elements very small, and then 
proceed to the limit. It was shown in Chapter 4 that the 
input current h 0 (t) of the artificial telegraph line is 

2 2< / it \ 

A o(0"£ c vrI °\crJ 

When the number of sections is made very large and the 
(' and R elements very small, the product CR becomes an 
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increasingly small quantity, and for values of * > o, the 
argument ( zt/CR ) of the /„ Bessel function becomes 
increasingly large. Thus we may *use the approximation, 


I 271 {cr]> f° r the / o 


Bessel function. 


Then the input current for the artificial cable reduces to 


A.W = 


J 


C _i 

R Vnt 


Equating this to the operational expression for the input 
current for the corresponding smooth line, we have, when 

C is made numerically equal to R , Vp^~y=, as before. 

Vnt 

The importance of these fractional order operators in 
transmission line theory is discussed in Chapter 7. 


6 



CHAPTER 7 

TRANSMISSION LINES 


TRANSFORMATION OF AN INTEGRAL EQUATION 
Consider a uniform transmission line having a distributed 
resistance R , a series inductance L, and a shunt capacity C, 
each per mile of length. For the moment we will neglect 
the effect of the distributed leakance G. In order to 
compute the characteristics of the travelling waves on such 
a system we must first formulate the indicial admittance 
function h(t). 

From alternating current steady-state theory we have 



where i represents the current x miles from the sending 
end of a uniform transmission line, of propagation con- 
stant P, terminated with its characteristic impedance Z 0 . 
In this expression V 0 denotes the sending-end voltage of 
frequency (o/m, and 


P- V jcoC(R + jtoL) 

7 I R +jwL 

°v jcoC 

In order to formulate the operational equation for h(t) } 
put V 0 « [i] and ja) - p. Then we obtain 

m. 

w \JR+pL 

This operational equation may be written in the more 
convenient form 

h(t) - [C. —JL = .g~f ' , J r +W 
\ L V> 2 + 2 <j)p 

where v - i/V LC, and <j> « R/zL 
74 



TRANSMISSION LINES 75 


Replacing the differential operator p by the parameter k, 
we may write the Carson integral equation as 


where 


" VF+ 2fX 




\ 



(II) 


The problem is thus reduced to evaluating the function 
h 0 (t) from the integral equation. The first step in finding 
a solution is to search through a table of Carson integrals 
for the required solution. If the required solution cannot 
be found in the tables it will be necessary to pick a solution 
that resembles the one sought and then transform it into 
the required form. Searching through the short table 
given in Chapter 6 , we find that the required solution is 
not given ; but the integral nearest the required form is 
found to be 

00 

U-*> V »+ 1 f 

VrfVl*-**)* ■ . (12) 


where J Q is a Bessel function of the first kind and t > b. 
In the integral the lower limit of zero may be replaced 
by b , since the Bessel function is zero for t<b, and exists 
only when t > b. It can be seen that the integral (12) 
resembles the integral equation (11) to some extent, and 
this resemblance suggests that we transform (12) in such 
a way that it will yield the function h 0 (t). 

Begin this transformation by putting Z = bk and t - bt x 
in equation (12). Making these substitutions, we obtain 
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In order to preserve our original notation in the parameter l 
we replace Z by X and t x by t. 

Now put A - o) + A, and obtain 
00 


e - ^(o>+h)'+b~ 

V(w+ij*+ 


" [*■ 
42 1 




Replacing co by A, and then making the substitutions 


A - ^ and 


we obtain 


e~V (■■ «+*)'+»* 




X 

V 


which may be written as 
^(tt+AiPW 





x 

r 




V , . 

where /q * -A and *» -J). In order to keep to our original 

notation we replace q by A and t % by t. Thus we obtain 
the Carson integral solution 






dt 


Comparing this solution with equation (11), it can be seen 
that they can be made identical by putting and 

b x -j(j>. Consequently the required solution may now be 
written as 


k_5Vj. +i# i 
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from which it can be seen that the function h 0 (t) of 
equation (n) is 

In this equation I 0 (x) = J 0 {jx), and is a Bessel function of 
the first kind with imaginary argument and of order zero. 
Consequently it follows that the indicial admittance func- 
tion h(t) and the solution of the operational equation may 
be written as 


m 


‘Jr 


Vp 2 + 2(j)p 

where v = i/VLC and (f> = R/iL. 


* + 2<f>p 




(i4) 


CHARACTERISTICS OF THE WAVE ‘ HEAD ’ 

Consider a point x miles from the sending end of the 
transmission line. From equation (14) it can be seen that 
h(t) is zero until t = x/v, at which instant it jumps suddenly 
to the value //, where 



After this initial jump h(t) may attain a maximum value 
which is large compared with the * head ’ H of the travelling 
wave as shown in Fig. 9. 

Rx 

As the length of the line x and the value of - — 

increases, the relative magnitude of the ‘ tail * as compared 
with the ‘ head ’ of the wave increases. When the length 
of the line becomes very great, or the value of L is decreased 
and becomes too small to be effective, the head H of the 
wave becomes negligibly small, and the wave, except in the 
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neighbourhood of its minute head, corresponds to the 
wave in the non-inductive ( L = o) line as shown in Fig. 9. 

The exponential factor in the equation for H is the most 
important, and in order to make the head of the wave as 
great as possible, the inductive time constant L/R should 



be made as great as possible in comparison with the capacity 
time constant CRx 2 . It can be seen that the exponential 
term is directly proportional to R and inversely proportional 
to Vl, hence it follows that a change in R makes itself 
felt in much greater proportion than one in L. 


THE EFFECT OF LEAKANCE 

Equation (14) neglects the effect of the distributed 
leakance G per mile. In order to obtain the operational 
equation for the indicial admittance when the distributed 
leakance G is included, replace pC by pC + G in the 
operational equation for h(t), and obtain 



-xP 
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where P « V(R + pL)(G + pC) 

The propagation constant P can be transformed into the 
form • 

, , R G 

where *'‘Tl + 2C 

, R G 

^“STaC 

and v denotes i/VLC, as before. 

The operational equation may be rewritten in the form 

- _ _J_( C + 

"vz&v + ?J' VifTji yzjf 

o 

which reduces the problem to an evaluation of the func- 
tion F(t). To evaluate this function we make use of 
the Carson integral solution expressed by equation (13). 
By putting h x » fa and b x - jfa we see that the explicit 
solution for h x (t) can be written as 


t 


‘‘-p)*- ■ fa) 


for 
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This equation is the exact solution for the current wave 
in terms of Bessel functions, and in order to evaluate the 
definite integral it is .necessary to expand the Bessel 
function as a power series, integrate term by term, or else 
evaluate it graphically. Thus it is a matter of some labour 
to calculate the effect of the leakance G irhany particular case. 

direct expansion of operational equation 
The foregoing solutions depend upon the transformation 
of the known Carson integral (12) into the required form : 
of equation (13). When such identities are known their 
value in connexion with the solution of the operational 
equations which arise in electric circuit theory requires no 
emphasis. Of course, we should not expect to find integral 
solutions suitable for use with every operational equation 
that is likely to arise. In such cases we can always fall 
back on the direct Heaviside method of expanding opera- 
tional equations. As an example of the direct method of 
expansion we will show how the function F(t ) can be 
deduced without the aid of Carson’s integral equations (11) 
and (13). 

From the formulation of h t (t) it can be seen that our 
problem is to evaluate the operator 

Shifting the exponential e~* li by means of the transforma- 
tion discussed in Chapter 6, we have 

m-'-** - p — 

p-K ^(p-h+w-p 

m -5 

This operational equation represents a function which 
x 

is zero until the time Nothing can happen at the 
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point x miles from the sending end until the wave arrives 
at x after a time interval t » x/v. This suggests that we 

transfer the time axis by means of the transfer operator e v . 
Making this transfer, the theory of which is discussed in 
Chapter 5, we get 


m 


-V* 

s Vtf'g y t g v 


Vpi-A ^p-iL.- -A 1 . 
P P P 2 p 2 a 2 lp a 


Vp 2 -^ 

We begin the development by expanding the term 
V/> 2 -<£ a a by the binomial theorem; we have 

3^ 6 

2 3 3 \p 5 

Let a denote x/v, and then form the infinite product 

*/ nt ,, , a*.* 3o^ t « 

a{p P *>)^ e 2p e 2'(2!)p* g 2'(3!)p“ 

Each of these expansions must be multiplied together and 
we then collect successive terms of inverse powers of p. 
Replacing p -n by t n /n\ in the operational expansion thus 
formed, we obtain 
*{p- t _ afy'ft t a 2 (j> 2 H 2 


=/(«) 

Thus we may write 

F(t)~ e-W.e-P a . 


(^h‘ + A 4 V+ 


-,-m 


_____ 

The next step is to expand the operator p/Vp 2 -^ in 
inverse powers of p, and substituting in the above expression 
for F(t), we obtain 

Fit) - j/(o + 


'• 3 v 5 /M“ 
2-4 ■f\p) j 


■ ■ • 


s 
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The next step is to perform the indicated integrations and 
gather the resulting series together in rising powers of t. 
After this the transfer operator e~ pa can be applied. 
Performing these operations, we obtain 



This result is due to Heaviside, who remarked that it was 
unfortunate that the series should be so complicated. By 
collecting terms in increasing powers of </> 2 and replacing 
a by x/v, Heaviside obtained the series solution 


m-c- 



which he recognized as the power series expansion of the 


function e~^ t l [ 


(s/-4 


Thus the function F(t) 


may be written from direct expansion as 


F(t)~e~m 


if 



which agrees with the result established by the Carson 
integral equation method. This example is an illustration 
of the fact that the direct Heaviside method of expansion 
may be tiresome to carry through. It must be noted, 
however, that only elementary mathematics are involved. 


Heaviside’s distortionless line 
From equation (15) it can be seen that the effect of the 
unit voltage applied at the sending end does not reach the 
point x of the transmission line until a time t - x/v has 
elapsed. Consequently v - x/t is the velocity of the wave 
and is due to the distributed inductance L and capacity C 
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of the line. For if we put R=G= 0, then (j> x - ^ - 0, 
and F(t) = e°I Q ( 0) - 1, and the indicial admittance func- 
tion h^t) reduces to 

x 

hj[t) = 0 for time t < - 

- rl for time t > - 
* 

which shows that the current wave jumps immediately to 
its steady-state value at time t - If an arbitrary 
electromotive force £(*) is applied at the sending end at 
reference time t = 0, then the resulting current wave i(t) is 
x 

i(t ) = 0 for time t < - 

- I-E(t- for time t > - 

VL \ » 

Consequently a non-dissipative transmission line transmits 
current waves without attenuation or distortion and with 
a finite velocity v. This non-dissipative transmission line 
is purely theoretical and quite unrealizable in practice. 
A distortionless line, however, may be realized physically. 
From a study of equation (15) Heaviside showed that the 
distortionless line is one in which 

, R G 
2L~2C =0 


and that 


x 

h(t) =» 0 for time t < - 


where 


~ rlg-fa f or time t >- 

V L * 

m * fC G II 
2\l L + 2 \/ C 


This shows that the current wave is transmitted without 
distortion of wave form. If an electromotive-force E(t) is 
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applied at the sending end, then the resulting current 
wave i(t) is 

x 

i(t) m 0 for time t < - 

" Jl e fXE ({~ f ° r time 1 > S 

These equations show that the distortionless line transmits 
the current waves without distortion of wave form, but 
attenuates the waves by the factor e ~ Thus the condition 

** 0 preserves the wave shape but introduces serious 
attenuation losses. It was one of Heaviside’s great achieve- 
ments to deduce the properties of the distortionless trans- 
mission line and to point out its approximately realizable 
character and to base thereon an engineering theory of 
telephonic transmission. 

Heaviside’s divergent expansions 

The series obtained for results so far have all been 
convergent. It was shown in Chapter 4 that these series 
are convenient for calculation only when the time t is small. 
When t is large, many terms need to be taken into account. 
Heaviside found that it is often possible to obtain easily an 
alternative form of series, which, in spite of the fact that 
it is divergent, can be used for computation when t is large. 
In order to fix ideas, let us consider an example; in 
Chapter 4 we discussed the Bessel function e'^I^t) and 
found that it was a simple matter to obtain a convergent 
power series expansion in the form 

«- . 1 - m + ^tf - + . . . 

This convergent solution is convenient for calculation when 
t is small. When t is large, however, this solution is 
difficult to use owing to its slow convergence. In order to 
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overcome this difficulty it is possible to obtain a divergent 
and asymptotic expansion of the form 




1 f i i 2 -3 2 I s . 3 s . 5 s 

Vtoft{ + (ip + 2i(8?i)‘ + " WW + ' ’ ' 


Now for relatively large values of t the successive terms at 
first become smaller, and as we proceed along the series 
we reach a smallest term. After this, however, the terms 
begin to increase and ultimately become very large. These 
characteristics are illustrated in Fig. 10. 



Fig. 10. — Graph showing effect of adding terms. 


Divergent series of the type illustrated in Fig. 10 are 
often convenient for calculation purposes provided we 
stop at the smallest term. Suppose, for example, that 
(8 <j>t) - 100, the second term in the brackets is O’Oi, while 
the third term is only 0-00045. The larger the value of 
(8 <f>t) the fewer the terms that need be used. For example, 
if (8 <j>t) > 100 the series reduces to i/V 2 n^t. In using 
such series for calculating the results of electric circuit 
problems we note that we must stop at the smallest term, 
beyond which it is useless to proceed. 

In order to show how these series arise in electric circuit 
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theory, put x = o in the operational expression for h(t) 
and obtain 


AoW' 


, ljc_ 

V R + p 


for the current entering the transmission line in response to 
a one-volt battery. There are several ways of solving this 
operational equation. Perhaps the quickest way is to search 
for the required form in a table of Carson integral solutions. 

Writing the expression for h 0 (t) as 


Jr-p 


‘ + PL 

note that the solution (E) in the short table of Carson 
integrals given in Chapter 6 gives the required form, 
namely, 


J 1 


4 




provided we put 2</> « R/L. From this result we see that 
the required solution is 




-pm ■ 


(16) 


where is a tabulated Bessel function. 

The equation (16) may be readily established by other 
methods. We may write 






2(f) 


where <j> - R/zL, and expanding by the binomial theorem, 
we obtain 

L3J/ 


Ht)- 


yi- H 


.LI 
2! \P. 


P + 


In order to convert this operational expansion into an 
explicit power series solution, we replace p~ n by t n /n\ in 
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accordance with the principles discussed in Chapter 3. 
The explicit solution is then 

*•»>- • ••} 

which is a convergent series solution of the operational 
equation in rising powers of /. This series solution is 
convenient for computing h 0 (t) when t is small, i.e. for 
calculating the first transient surge in electric circuit 
problems. When t is large, however, this series is difficult 
to use owing to its slow convergence. In order to overcome 
this difficulty it is necessary to recognize and sum the series. 
After a trial it will be found that the series is the convergent 
power series expansion of the tabulated Bessel function in 
equation (16). If we are unable to recognize and sum the 
power series expansion and wish to calculate h 9 (t) for large 
values of t , it will be necessary to use the Heaviside scheme 
of expanding the operator in rising powers of p , from 
which a convenient series for calculating purposes can be 
obtained. The operational equation for h 0 (t) may be 
written in the form 



Expanding the denominator by the binomial theorem, we 
get 



It will be observed that the operand is now V/>[i], and 
consequently we make use of the fact, rigorously established 
in Chapter 6, that V^fi] must be interpreted as 1 /\/nt. 
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d n 


Making this substitution and replacing p n by we obtain 


h 0 (t) = 



i d 1.3 1 d 2 

rf'dt * IT '(4^'A 2 


1.3.5 i d 8 | 1 1 

~T w^ + ■ ■ 


Performing the indicated differentiations, we obtain 

• 3 - 5 )\ 


'«(*) 


V L Vzmfrt \ 


(1.3 Y (1 


This solution may be used to compute h 0 (t) when t is large, 
and is, in fact, the asymptotic expansion of the function 
e ~ The asymptotic expansion of this function is 

usually derived by intricate and difficult mathematical 
processes and the remarkable ease of the Heaviside method 
should be noted. 


cable with terminal resistance 
As an example of Heaviside’s method of obtaining 
two expansions (one for small values of t and the other for 
large values of t), consider the problem of the determination 
of the voltage e 0 (t) which appears at the terminals of a non- 
inductive cable when a one-volt battery is applied through 
a terminal resistance R 0 . The arrangement is shown in 
Fig. 11. 



R ./2 


Fig. 11.— Non-inductive cable with terminal resistance. 
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The impedance function of the resistance R 0 and the 
cable in series, is 

where R = resistance per mile of cable, 

C - capacity per mile of cable. 

The entering current h 0 (t) is given by 


/>»(/)= 


R a + V R/pC 

and the voltage e 0 (t) at the cable terminals is 
Ra 


- 


R 0 + VR/pC 


(17) 


(18) 


In order to obtain a convergent power series solution we 
must expand the operational equation (17) in inverse powers 
of p. Writing equation (17) in the form 


h 0 ( 0 - 


R 0 'i+ i/(hp)t 


where 


Vh = R 


7 : 


We obtain on expanding and collecting terms 
1 


K(t)- 


1+ 


«.r m w 


(MJ- i + _l + 

■' pw ‘ ’ ‘ 


Replacing p -B by i™/n! and Vp by i / Vjtf, we obtain 






_\^j 1 , 1 1 

~R„V*l(W + (W + --j 

Performing the indicated integrations in the second series 
7 
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and recognizing the first series as the expansion of the 
exponential we obtain 



Substituting in equation (18) we obtain for the terminal 
voltage e 0 ( t ) the expression 


«»(<) 



This result is convergent and convenient for calculation 
when t is small. Although this series is valid for any value 
of t it is hard to compute from when t is large. It will be 
noted that the voltage e Q ( l ) contains a positive exponential 
term ; this voltage, however, does not increase expon- 
entially, because for large values of f, the exponential term 
is practically cancelled by the series that follows it. In 
order to obtain an asymptotic expansion suitable for calcula- 
tion when t is large, we must expand the operational 
equation (17) in ascending powers of p. From equation (17) 
we have, after a simple rearrangement 




(W 


R 0 ‘i + (hp)t 
We obtain on expansion 
1 


where Vh~ R { 


■ji 


K(t) - F (hp)i{i - (hp)i + (hp) - (hpy + (hpf - 


Now for the finite values of the time we are concerned ' 
with, p n [i] - 0, where n is any positive integer. Thus the 
integral powers of p may be dropped, and rearranging we 
have 
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Replacing Vp by i/y/nt, and performing the indicated 
differentiations, we obtain 



Thus, the voltage e 0 (t) across the cable terminals is 



which is an asymptotic expansion convenient for com- 
putation when t is large. It must be noted that the error 
involved in using an asymptotic series does not become 
negligible until / becomes large. This, of course, is all 
that can be expected from an asymptotic series ; if t is 
small it is a simple matter to use the ordinary convergent 
series. 

This example is an illustration of Heaviside’s method of 
obtaining divergent and asymptotic series solutions of 
electric circuit problems. Heaviside found that if an 
operational equation, h(t) - i/Z(p\ admits of series 
expansion in the form 

h(t ) m Co + CjVj p + C 2 p + CipVp + C t p* + . . . 


a solution, usually divergent and asymptotic, results from 
discarding the terms in integral powers of p , and replacing 
d n 1 

p n Vp by Thus the explicit solution may be 

written 

W-c.+jc.+c^+c.**..]^ 


where the C’s are constants. 

It can be seen that Heaviside’s asymptotic expansions 
may be of great value , in certain kinds of electric circuit 
investigations. The question of when and how divergent 
and asymptotic series may be rigorously developed involves 
advanced mathematical considerations. Heaviside did not 
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attempt to explain the mathematical theory; the only 
justification he gave was an appeal to direct measurement 
to check his calculations. He appeared to have evolved 
a simple working system based on trial, confirmed on test, 
and then applied to electric circuit problems on a large 
scale. By virtue of his last theorem, however, it would 
have been an easy matter for Heaviside to have investigated 
the mathematical validity of his operational methods. 
Thus if he split his operator F(p) into a sum of terms and 
if when he evaluated each term by his integral theorem they 
added up to F(A)/A, then the operational process would be 
justified : there are, of course, other methods. It is 
possible that Heaviside made these investigations but never 
bothered to publish them ; in this case they will probably 
be found reposing peacefully in the archives of the Institu- 
tion of Electrical Engineers together with other unpublished 
work of Heaviside. Since the time of Heaviside, however, 
the theory of his divergent and asymptotic expansions has 
been placed on a rigorous mathematical footing by Wiener 
and Carson. 



CHAPTER 8 


THE APPLICATION OF MODERN THEORIES OF 
INTEGRATION TO THE SOLUTION OF 
CIRCUIT PROBLEMS 

COMPLEX INTEGRATION 

In Chapter 5 it was stated that the Heaviside operational 
equation is not a true algebraic equation at all but is simply 
a shorthand way of writing down the descriptive differential 
equations of the netwprk. It was also shown that to every 
operational equation there corresponds a Carson integral 
equation. This integral equation expresses a unique but 
implicit relationship between h(t) and Z(X), which enables 
us to extend the range of operational solutions by direct 
transformations. The formal solution of Carson’s integral 
equation is given by the Bromwich contour integral which 
gives an explicit relationship between h(t) and Z(A). By 
virtue of the Carson-Bromwich relationships it is possible 
to establish rigorously all the fundamental theorems of the 
Heaviside operational calculus. 

It may seem curious to the student that the operational 
calculus, invented by Heaviside for dealing with functions 
of a real variable, should find it necessary to use the theory 
of functions of a complex variable in order to establish 
rigorously its fundamental theorems. This is because the 
descriptive differential equations which arise in network 
problems are expressible in terms of analytic functions. 
These analytic functions are found to satisfy Cauchy’s 
integral theorems. Consequently use is made of the 
calculus of residues in order to evaluate the infinite integrals 
which arise. 

In the theory of functions of a real variable we consider 
functions associated with the movement of points along a 
straight line, and begin the study of these functions by 

93 
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considering the class of functions which possess derivatives 
of all orders and are capable of expansion in a power series 
by Taylor’s theorem abqut any point on the line. In the 
theory of functions of a complex variable, however, we are 
dealing with functions associated with movements of points 
in a plane, and begin the study of these functions by con- 
sidering the restricted class of analytic functions which 
possess derivatives of all orders and are capable of expansion 
in a power series by Taylor’s theorem about any point in 
a given region in the plane. Thus if f(z) is analytic in 
a given region in the complex plane it can be expanded in 
a Taylor series about any point z - a in that region and 

f(z) = Za n {z-a) n 
0 

In order to apply complex variable theory to electric 
circuit problems we begin by considering the meaning of 
the integral of a function of a complex variable along a plane 
curve. Now the physical ideas associated with line integrals 
in the complex plane are similar to those associated with 
integrals taken along the real axis, except that the integrand 
is of the form x+jy } and consequently the differential is 

also complex. For example, the meaning of |/(^) dz 

between any two points in the complex plane may be 
visualized as follows : the path of integration may be 
divided into tiny segments dz , and this elemental length is 
multiplied by the value of the function/^) at the midpoint 
of the segment. We then add up all these elemental 
products and take the limit of the sum as <5* o. This 

limiting sum is then the value of the integral 1 f(z)dz 

between the two points considered. It follows from these 
considerations that the integral around a closed path is 
zero, provided there are no singularities or points of dis- 
continuity on or within the path of integration. This 
statement concerning the value of f(z) taken round a closed 
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contour is the basis of a fundamental theorem due to 
Cauchy, which may be stated formally as follows : if 
f(z) is analytic at all points on and inside a closed contour C, 

then J* f(z) dz = o. Because of the fact that the line integral 

taken a closed contour C is zero, it follows that the path 
of integration may be stretched or contracted in any way 
we please, provided in so doing we do not pass over any 
singularities of the function, that is, we must not pass over 
points at which f(z) becomes infinite. Should the area 
contained within the path of integration enclose a singular- 
ity, then we may contract the path until parts coincide and 
cancel, and we are left with only a small circle surrounding 
the singularity. These ideas are illustrated in Fig. 12, 
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where C is a contour which contains a singularity or point 
of discontinuity a. 

It can be seen from Fig. 12 that if the contour C is 
contracted so that the parts coincide and cancel, we are 
left with only a small circular contour surrounding the 
singularity at a. 

These ideas are important in electric circuit theory ; for 
example, in the examination of a circuit problem by means 
of Heaviside’s expansion theorem we see that the roots 
of the equation, Z(X) = 0, play an important part. For 
dissipative systems these roots are complex with negative 
real parts, and, if plotted in the complex plane, they will 
generally lie on the left-hand side of the axis of imaginaries. 
If, however, the system contains negative resistances, then 
it is possible that the complex roots may have positive real 
parts. These roots are related to the singularities of the 
function i/Z(A) and determine the nature of the operator 
with which we have to deal and the process which may be 
safely employed to interpret it. 

If/(;r) is analytic in a given domain in the complex plane, 
except at the singularity 2= a, then we can draw two 
concentric circles of centre a , both lying within the domain. 
In the annulus between these two circles, we have by virtue 
of Laurent’s theorem, an expansion of the form 

p) = Ea n (z - af + Zb n (z - a) 

0 1 

The second term on the right-hand side is called the 
‘ principal part ’ of f(z) at the singularity z * a. It may 
happen that b n + 1 = b n+2 = £ n+3 - . . . = 0 but b n + 0. 
In this case the principal part consists of the finite number 
of terms 

»i , *» , . 

(z-a) (z-af ' ' ' (z-aY 

and the singularity at 2 - a is called ‘ a pole of order n ’ 
of fix) and the coefficient b x of i/(z- a) is called ‘the 
residue of fa) at the pole z - a \ If the principal part is 
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an infinite series the singularity at z = a is called ‘ an 
essential singularity \ In electric circuit theory the func- 
tional singularities are usually first-order poles or essential 
singularities. From the above definition it can be seen 
that b u the residue of the function at a first-order pole, is 
given by b t = lim (z - a)f(z). 

z-y a 

It follows from Cauchy’s theorem that the value of a line 
integral about a closed contour is equal to the sum of the 
integrals about small circles in the original area each 
enclosing a pole. It is an easy matter to show that the line 
integral of a function taken around one of these small circles 
is equal to 2 nj times the residue b l of the function at the 
pole. For example, suppose we integrate the simple 
function i/A about one of these circles. It will be noted 
that this function has a pole at the origin and consequently 
the line integral of i/X about the origin is 21 rj } for the 

residue there is limA[^ j which is unity. If, however, 

W 

the circle did not enclose the pole at the origin, the value 
of the line integral of i/X would be zero. This can be 
easily seen by integrating round the circle shown in Fig. 13. 



Fig. 13, — Integration of line integral round circle. 
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Putting l • pe jd , we have 

|j = PogA] c -[Iog | A | +jO]c 

c 

(change in 0) 
c 

It can be seen that the change in 6 is 0 if r < a and 
in if r > a> which is in agreement with the result established 
by the residue method. The statement that the value of 
a line integral round a small circle enclosing a pole, is 
equal to 27 tj times the residue b u of the function at the 
pole, is a special case of Cauchy’s Residue Theorem which 
may be stated formally as follows : if f(z) be analytic on 
and within a closed contour -C, except for a finite number 
of poles within C, then 

j /(*) dz — 2 njER 

c 

where ZR is the sum of the residues of f(z) at its poles 
within C. In the majority of electric circuit problems the 
poles which arise are mostly first-order poles and conse- 
quently the residue at a pole of f(z) may be calculated 
from lim (z - a)f(z). 

z-ya 

The application of Cauchy’s residue theorem to electric 
circuit theory is important because it will enable a student 
to see clearly the functional aspects of Heaviside’s opera- 
tional calculus. It will enable him to extend the method 
to unusual problems and to see danger signals in doubtful 
cases. He will also find that an argument in terms of 
integrals in the complex plane is very much simpler than 
one in terms of ordinary integrals. 

THE BROMWICH CONTOUR INTEGRAL 

We have yet to demonstrate that the Bromwich contour 
integral (9) is the formal solution of the Carson integral 
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equation (8). There are several ways bf doing tliis ; we 
recall that the function h(t) is uniquely but implicitly 
determined by the Carson relationship, and consequently 
the problem is to formulate h(t) as an explicit function of 
Z(X). Since we are primarily interested in the application 
of operators to a discontinuous operand, we begin by 
formulating the unit function [i] as a contour integral. 
This unit function is zero before, unity after time t = o, 
and is consequently discontinuous at the origin. From 
Cauchy’s residue theorem it is apparent that the unit 
function can be expressed as 

(’ 9 ) 

C 

provided we choose a suitable contour and a function 
of X that is everywhere regular in the X plane, except at 
the origin, where it must be discontinuous and have a 
residue of unity. It is a simple matter to see that such 
a function is e n /X. This function has a pole at the origin 

f e u\ 

in the X plane and its residue there is lim X [ — ), which 

is unity. Consequently f(X) in equation (19) may be 
taken as e H /X, and C may be a circle surrounding the 
origin in the X plane as shown in Fig. 14. 

We know that it is the reciprocal of the impedance 
function that acts on the unit function and that the complex 
roots of Z(X) = 0 have negative real parts for all dissipative 
systems. We must therefore see that all the singularities 
of the integrand are on the left-hand side of our contour. 
Consequently we enclose the origin in the X plane by a 
small half-circle so as to include it in the left-hand semi- 
circle as shown in the upper part of Fig. 14. 

Now the value of the line integral taken round the right- 
hand contour ABBA in the upper part of Fig. 14 is 
obviously zero because the path of integration is regular 
and the enclosed area contains no poles. The value of 
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the line integral taken round the left-hand contour ADEA, 
however, is z nj times the residue at the origin, which in 
this case is unity. Consequently it is possible to write 
equation (19) in the form 

<”> 

V 

with the contour shown in Fig. 14. An equivalent path 



would be a line parallel to the imaginary axis, but at a 
constant distance C to the right as shown in Fig. 14. 
These considerations led Bromwich to write the Heaviside 
operational equation (4) in terms of the contour integral 

C+joo 



C-j 00 


where C is so chosen that all the singularities of the 
integrand lie on the left of the path of integration in the 
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complex plane. In the majority of electric circuit problems 
the constant C may be taken as zero ; there are exceptions, 
however, and cases may arise when the Bromwich equa- 
tion (9) is valid only when C is made greater than some 
finite constant. Such an exceptional case arises when the 
electrical system specified by Z(X) is ‘ unstable ’, that is, 
it may contain some internal source of energy or ‘ negative * 
resistance which makes its transient oscillations increase 
with time. In such a case, Z(X) = 0 will have roots to 
the right of the imaginary axis, and in order that the 
Bromwich contour integral shall always give the correct 
interpretation of its corresponding operational equation, 
C must be large enough to ensure that all the singularities 
of the integrand lie on the left of the path of integration. 

DERIVATION OF EXPANSION THEOREMS 

In Chapter 2 we showed that the function i/Z(X) can 
be expressed in terms of partial fractions, and that 

1 y > .. 

Z(2J Z(A- WJg) 

Substituting this result in the Bromwich equation (9), we 
obtain 

c-tioo n 

Kt)m k J ■ (2,) 

c-i 00 

In order to apply Cauchy’s residue theorem to the evalua- 
tion of this integral, we first note that the integrand has 
poles at 0, X ly A 2 , . . . X K , . . . X n . Some of these poles 
and the surrounding contours are shown graphically in 
Fig. 14. The contour on the left-hand side includes 
within its area the poles at 0, X ly X it . . . X IL . We may now 
contract the contour ABCDA in such a way that we are 
left with only small circles surrounding the poles as shown 
on the right-hand side of Fig. 14. We must now evaluate 
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the residues at the poles of the integrand in equation (21). 
The residues at the pole at the origin of the integrand is 

which is the steady-state term. 

The residue at the pole 1 K of the integrand is 

Rk 

i5k) 

Hence the sum of the residues is 

n 1 i 
_ x ri e K K l 

ZR -z(o) + Zr K zW) 

and by virtue of the Bromwich formulation we may write 

n 1 t 

u \ 1 V 7 e * Rt 

which is the well-known Heaviside expansion theorem. 
This theorem has now been deduced by a rigorous mathe- 
matical method and verifies the simple experimental method 
used in Chapter 2. 

The partial fraction expansion used in this derivation 
shows that the roots of the impedance function, Z(A) = 0, 
must be unequal and distinct from each other. In all 
finite dissipative networks this condition is satisfied. The 
case of repeated roots can be dealt with by the method 
used in Chapter 2, namely, letting the roots approach 
equality as a limit. 

Tins method of derivation may be easily extended to 
cover alternating voltages ; suppose, for example, that 
instead of the unit voltage [1], we applied the alternating 
voltage E cos cof[i] to the network at time t » 0. Let it 
be required to determine the current response i(t) of a 
network described by the impedance function Z(X). For 
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purposes of analysis we may write this alternating voltage 
as the real part of Ee i<at [ 1] ; and from the Bromwich 
equation (20) it can be seen that this type of voltage may 
be represented by the contour integral 


EeW[ 1] * 




™j\x 


M 


with the same contour as for equation (20). It follows 
from our discussion of the differential equations of circuit 
theory that if the applied voltage is of the form, 
then Z(k+j(o) will represent the impedance function of 
the network. Thus we may write i(t) as the Bromwich 
contour integral 

G+j 00 

e (nju)t 

: dl . . . (22) 


m 


E_ 

2 Jlj 


XZ(X+jco) 


C-jo 0 


The integrand has poles at 0, (A x - /ft)), (A 2 -/w), . . . 
“>), • • • 

The residue R 0 at the origin is 

D , Ee^+W EeW 

which is the steady-state term. 

The residue R K at the pole (h K -jo) is 

• Ee^+W 

Rk ' (A * ~ ] 0 >)} -X{X-(X k -W\Z\X^J>) 

Eeh 1 

RK ~{X K -jm)Z'{X K ) {K ~ ■’ 2 ' 3 ’ • ■ • n) 
Hence i(t), which is the sum of the residues, may be written 


Ee 1 w t 

* )m Z(M + E 2 ^ i (X K -j<n)Z'(X K ) 

which is the alternating-current expansion theorem dis- 
cussed in Chapter 3. The first term on the right-hand 
side of the alternating-current expansion theorem represents 
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the steady-state current, while the summation terms give 
the transient oscillations. The alternating-current expan- 
sion theorem was given by Carson ; it will be noted that 
when ca » o and E - i, Carson’s expansion theorem 
reduces to Heaviside’s expansion theorem as it should do. 

A RELAY PROBLEM 

The Carson expansion theorem is sometimes very useful 
for the solution of alternating-current transient problems 
which arise in connexion with the design of lumped 
circuits. As an example, suppose it is required to investi- 
gate the transient oscillations in the circuit shown in Fig. 15. 



At reference time t = 0, a direct voltage E B and an 
alternating voltage E cos c ot are applied to a circuit con- 
taining resistive and inductive elements as well as a charged 
condenser. It is required to investigate the transient 
oscillations in this circuit, which is part of a relay arrange- 
ment used in automatic telephony. 
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From the principle of superposition we know that at 
time t > 0, the current i(t) in the circuit will be 

i(t) - h(t) + H{t)\ h(t) 

where i x (t) = current due to discharge of condenser 
h(t) - current due to direct voltage E B 
i 3 (t) = current due to alternating voltage. 

In order to obtain an expression for i x (t) we note that if 
a battery of voltage E 0 be suddenly connected in series 
with an uncharged condenser, the resulting current is the 
same as when the condenser is initially charged to a 
voltage E 0 and then discharged through the circuit, the 
direction of the voltage being maintained the same in the 
two cases. This follows because in each case the voltage e(t) 
impressed on the circuit is initially E 0 and is decreased by 

an amount, gji dt, as the current flows. Since the voltage 

impressed on the circuit is the same in each instant, the 
same current will flow, and one may be analysed in place 
of the other. Consequently the problem of an initially 
charged condenser is reduced to that of a direct voltage 
applied to an uncharged condenser. Thus the operational 
equation for i\(t) may be written as 

R + pL+i/pC 

where i? « + R 2 + R z 

E 13 Li + 


Also E 0 may be replaced by Q/C> where Q is the initial 
charge on the condenser in coulombs and C is the capacity 
in farads. It will be recalled that the solution of this 
operational equation was obtained in Chapter 2. We are 

1 /? 2 

here interested in the case where — > — ; thus we 
may write 


8 


h {t)-^ c .e-«s m sh 
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where a - R/zL and ft » 

The operational equation for t 2 (<) is the same as for ^(t), 
except that E 0 is now replaced by E B , the direct battery 
voltage. Thus we may write 

h(t) = ~ e ~ oU sm^ t 

An expression 'for i z (t) can be obtained by applying the 
Carson expansion theorem. The values of A 2 and 
dZ(X)/dX are, of course, the same as those obtained in 
Chapter z for the impedance function Z(p) = R+pL+ i/pC. 
Substituting for i t (t) we obtain 

R+jojL+i jcoC zfiL j/o+a-j/ joj+ct+jp J 

Simplifying this expression and then taking the real part, 
we obtain 

i ~ 4) + ( a2 + ^ 2 ) s * n (0* * v) 

" VF-TX* - + 7^ ' 

where tan (f> = AT//? 

X = wL- i/oC 

tan y> - 2 ft) 2 ot/?/ {(a 2 + jft 2 ) 2 + co 2 (a 2 - /9 2 )} 

For a) = o, the steady-state component becomes zero and 
the alternating-current transient component reduces to the 
direct-current transient. If the transient distortion D be 
defined as the instantaneous difference between the actual 
response i(t) of the relay circuit and the steady-state 
response, we see that the transient distortion may be 
written as 


/_i_ if 2 

v LC~ 4L 2 



l(«» + fl 2 )8in (fr-y) ] 
V (a 2 + ft 2 - ft) 2 ) 2 + 4a 2 ct) 2 J 


Now 0/C is the initial voltage E 0 across the condenser at 
reference time t * 0 , and it can be seen that if the battery 
voltage E b is made equal and opposite to E 0 , the first term 
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on the right-hand side of this equation will be zero. This 
shows that the transient distortion in the relay circuit may 
be considerably reduced by making E B oppose E 0 . For low 
Frequencies the second term on the right-hand side of the 
equation is small. 

THE SUPERPOSITION INTEGRAL 
It will be observed that in the relay problem we have 
just discussed, i(t), the resultant current, was written as the 
sum of the component currents t\(t), i 2 (t) and i 3 (t), by the 
principle of superposition. This principle states that when 
several electromotive forces are applied simultaneously to 
a network of fixed parameters each produces its own effect 
independently of the others. Thus the result of several 
electromotive forces acting simultaneously is found by 
adding the effects which would be produced by each act- 
ing alone. The principle of superposition was used in 
Chapter 1, where it was stated that when the indicial 
admittance function h(t) has been formulated from the set 
of auxiliary differential equations (2), it is possible to 
obtain solutions of the canonical equations (1), in the form 

t 

d f 

'»(')=■ J «('-#»( v) 4 • • ■ (3) 

0 

where i n (t) represents the current which will flow in the 
nth mesh of the network in response to an impressed 
voltage e(t) of arbitrary wave form. Equation (3) can be 
established in a number of different ways. For example, 
the current i(t) may be written in operational form as 
i(t)-e(p)h(p) 

where h(p) is the operational equivalent of h(t). By virtue 
of Heaviside’s last theorem we may write h(p) as 

00 

h(p) * |a(v>)£ ~ m p dy) 

0 
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and consequently 

oo 

e(p)h(p) - p^e - plf e(p)h(rf) df 

0 

Replacing e(p) in the integrand by its equivalent time 
function and performing the necessary transfer operation, 
we may write 

t 

>W - 

0 


which is the same as equation (3). If we make the sub- 
stitution x ~ t-ip, and then in order to preserve' our 
original notation, we change the variable of integration 
from x to y), we may write equation (3) in the form 


t 



0 


It will be noted that this formula contains a variable 
parameter under the sign of integration, and in order to 
differentiate this expression between the variable limits t 
and 0, we put 

t 

0 

t 

then j,# : Y>) “ jj/K* ~ vKv) d f + h (°W) 

0 

t t 

Thus h(t - f)e(y}) d\p = h(o)e(t) + |a'(< - tp)e(y)) dy) 

0 0 

Likewise we have 




MODERN THEORIES OF INTEGRATION IO9 

where the primes denote differentiation with respect to the 
argument. This superposition formula enables us to find 
the response of a network to an applied electromotive force 
which varies in any manner whatever, provided we know 
the value of h(t). The first term on the right-hand side is 
the current term which results when the switch is closed, 
and is zero if the voltage passes through the origin at time 
t - 0. The effect of the remainder can be added to this 
by finding the effect of a step function which closely fits 
the function e(t ) and then passing to the limit. If the 
superposition formula in any particular case is difficult to 
evaluate formally, then it is always possible to evaluate it 
by graphical or mechanical processes. On this fact rests 
the practical as distinguished from the purely theoretical 
value of the formula. 


THE FOURIER INTEGRAL 

If it is found that the steady-state solution of an electric 
circuit problem is sufficient, then the need for the use of 
the Heaviside operational calculus discussed in this book 
does not arise. In cases such as this the well-known 
symbolic method is sufficient. Perhaps it is now generally 
recognized that it is not possible to predict completely the 
performance of an electric network from a single steady- 
state solution. It is necessary to integrate the individual 
steady-state responses over an infinite frequency range. 
Now it is important for the student to realize that there is 
no fundamental mathematical difference between the 
Heaviside ‘ operational ’ method of solution and the Fourier 
‘ steady-state ’ summation method. Both methods are ways 
of solving the canonical equations of electric circuit theory 
discussed in Chapter 1 . It can be shown that the Heaviside 
Calculus is a simple practical method of evaluating the 
infinite Fourier integrals which arise in electric circuit 
problems. 
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The Fourier solution of the canonical equations discussed 
in Chapter i may be written 

00 00 

if r 

J E{v)dv J -m M - ■ (23) 

- 00 - 00 

This infinite integral formulates the current response i(t) of 
the network, specified by the impedance function Z(jco), 
when an electromotive force E(t) is applied at time t - o. 
The general method of evaluating this Fourier integral is 
by means of contour integration in the complex plane and 
the calculus of residues. By this process it may be applied 
to the solution of electric circuit problems ; but, generally 
speaking, this method of solution is very difficult to apply. 
Compared with the Heaviside method, the Fourier integral 
solution (23), has no advantages from the standpoint of 
rigour, and presents formidable mathematical difficulties 
(even to professional mathematicians) when direct evalua- 
tion is attempted. 

It is important to note that the Carson integral theorem 
may be derived from the study of electrical networks in 
terms of the Fourier integral, and is, in fact, a general 
relationship resulting from a simple application of the 
theory of Fourier transforms. In this theory, if 

jeo 

-JOO 

00 

then <Hp o)= j e-«i (h)dX 

- 00 

is called the Fourier transform of f(t). Now suppose the 
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indicial admittance function h(t) be written as the Fourier 
integral 


j oo • 



-;co 


then the Fourier transform of this gives the Carson integral 
relationship 

00 

/(«) - <uJc-" A A(A) dk 

0 

since h(t) = o when t < o. 

Thus it follows that the operational methods associated 
with the Carson integral theorem are the working tools of 
the classical Fourier analysis. 

Although the direct evaluation of the Fourier integrals 
which arise in the analysis of electric circuit problems 
presents formidable mathematical difficulties, it must be 
noted that it is possible to use the Fourier integral as a 
means for obtaining prafctical design criteria without getting 
involved in the mathematical pitfalls associated with direct 
evaluation. For example, in an electrical communication 
system, the network must be so designed that the received 
signal is a more or less faithful copy of the sent signal. 
The departure from this requirement is a measure of the 
distortion introduced by the network. The Fourier integral 
provides a simple method of isolating the factors on which 
this distortion depends and of formulating design criteria 
for distortionless transmission through the network. Thus 
it can be shown without the explicit evaluation of the 
Fourier integrals concerned, the conditions for distortion- 
less transmission through the network are that over the 
essential range of frequencies contained in the sent signal, 
the transfer impedance of the network be equalized both 
as regards amplitude and phase ; that is, the amplitude 
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must be constant and the phase angle linear with the 
frequency. For another example it is interesting to con- 
sider the formulation of the indicial admittance function h(t) 
as a Fourier integral. From a study of this integral it can 
be shown that the behaviour of a network under all circum- 
stances is completely determined if either the real or 
imaginary component of the admittance function i/Z(;w) 
is specified over the entire frequency range. The direct 
evaluation of h(t) from the Fourier integral formulation is 
a very complicated procedure, and unless the infinite 
integral can be recognized, the only practical method of 
evaluating h(t) is by means of the operational calculus 
discussed in this book. 

In Chapter 5 Carson’s integral theorem was obtained as 
a corollary to Heaviside’s last theorem. It is a simple 
matter to show that by virtue of these theorems the rules 
and formulae of the operational calculus may be simply but 
rigorously deduced. Heaviside, however, did not use these 
theorems to establish operational formulae and processes ; 
but employed them to evaluate the infinite integrals which 
arise in electromagnetic problems. The theorems are of 
use either way. Since there are other valid ways of estab- 
lishing operational formulae, the use that Heaviside made 
of his last theorem is probably the most valuable. It is 
sometimes an exceedingly difficult job to evaluate an infinite 
integral without the aid of this theorem. In the third 
volume of his Electromagnetic Theory Heaviside gives some 
seventy examples of the evaluation of integrals by the 
application of the impulse theory discussed in Chapter 5. 
In these examples he discusses not only the evaluation of 
ordinary infinite integrals but shows how the fundamental 
Fourier theorems may be readily established by the 
application of his impulse ideas. He also deals with Bessel 
theorems, elliptic functions and, among other things, with 
convergent and divergent series. In all cases his opera- 
tional treatments are interpreted in accordance with the 
impulse ideas discussed in Chapter 5. As a simple example, 
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00 

suppose that the infinite integral xt) dx is en- 

o 

countered. Replacing the Bessel function J 0 (2V xt) by its 

_ X 

operational equivalent e />[i], and then integrating with 
respect to x> we obtain p[i]. The value of the infinite 
integral is therefore zero, except for the special case where 
t - o, in which case it is infinite. This result, which has 
been obtained so easily, would be very difficult to establish 
by other methods. Only this one example will be reviewed 
here, for our subject is Heaviside’s electric circuit theory 
and not the evaluation of infinite integrals by special 
methods invented by Heaviside. A complete discussion of 
the evaluation of infinite integrals by means of Heaviside’s 
impulse functions is beyond the scope of an elementary 
book such as this. 
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